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ABSTRACT 

We present the serendipitous discovery of molecular gas CO emission lines with the IRAM Plateau de Bure 
interferometer coincident with two luminous submillimeter galaxies (SMGs) in the Great Observatories Ori- 
gins Deep Survey North field (GOODS-N). The identification of the millimeter emission lines as CO[4-3] 
at z = 4.05 is based on the optical and near-IR photometric redshifts, radio-infrared photometric redshifts 
and KeckH-DEIMOS optical spectroscopy. These two galaxies include the brightest submillimeter source in 
the field (GN20; ^gso^m = 20.3mJy, zco = 4.055 ±0.001) and its companion (GN20.2; ^gso^m = 9.9mJy, 
^co = 4.051 ±0.003). These are among the most distant submillimeter-selected galaxies reliably identified 
through CO emission and also some of the most luminous known. GN20.2 has a possible additional counter- 
part and a luminous AGN inside its primary counterpart revealed in the radio. Continuum emission of 0.3mJy 
at 3.3mm (0.65mm in the rest frame) is detected at Sa for GN20, the first dust continuum detection in an SMG 
at such long wavelength, unveiling a spectral energy distribution that is similar to local ultra luminous infrared 
galaxies. In terms of CO to bolometric luminosities, stellar mass and star formation rates (SFRs), these newly 
discovered z > 4 SMGs are similar to z ^ 2 — 3 SMGs studied to date. These z 4 SMGs have much higher 
specific SFRs than typical B-band dropout Lyman break galaxies at the same redshift. The stellar mass-SFR 
correlation for normal galaxies does not seem to evolve much further, between z 2 and z 4. A significant 
z 4.05 spectroscopic redshift spike is observed in GOODS-N, and a strong spatial overdensity of B-band 
dropouts and IRAC selected z > 3.5 galaxies appears to be centered on the GN20 and GN20.2 galaxies. This 
suggests a proto-cluster structure with total mass ^ IO'^^Mq. Using photometry at mid-IR (24/im), submm 
(850/xm) and radio (20cm) wavelengths, we show that reliable photometric redshifts (Az/ (1 +z) 0.1) can 
be derived for SMGs over 1 ^ z ^ 4. This new photometric redshift technique has been used to provide a first 
estimate of the space density of 3.5 < z < 6 hyper-luminous starburst galaxies, and to show that they both 
contribute substantially to the SFR density at early epochs and that they can account for the presence of old 
galaxies at z 2 — 3. Many of these high-redshift starbursts will be within reach of Herschel. We find that the 
criterion 4 ghz ^ S2/i,^m^ coupled to optical, near-IR and mid-IR photometry, can be used to select z > 3.5 
starbursts, regardless of their submm/mm emission. 

Subject headings: galaxies: formation — cosmology: observations — infrared: galaxies — galaxies: star- 
bursts — galaxies: high-redshift — submillimeter — galaxies: individual (GN20, GN20.2a, 
GN20.2b) 



L INTRODUCTION 

Submillimeter and millimeter observations of galaxies are 
subject to strong negative K-corrections over a very broad 
range of redshifts, and can reach nearly constant bolometric 
luminosity limits (for a given flux density limit) at 1 < z < 7 
(Blain et al. 1993). This provides a powerful method for de- 
tecting and studying starburst galaxies in the distant universe. 
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Efforts to identify counterparts of bright (^gso^^m > 5 mJy) 
submm-selected galaxies (SMGs; see Blain et al. 2002 for a 
review) based on radio detections have recently suggested that 
the typical redshift of these sources is z ^ 2 (Chapman et al. 
2003; 2005), with few objects at z ^ 3 and virtually none at 
z > 4. This result suggests a very strong decline in the space 
density of the highest luminosity, most highly star forming 
galaxies, in the distant universe. 

SMGs are outliers in the stellar mass-star formation rate re- 
lation in the z ~ 2 universe (Daddi et al. 2007a; Takagi et 
al. 2008), have gas kinematic signatures suggesting ongoing 
mergers (Tacconi et al. 2006; 2008), have very compact sizes 
reminiscent of massive, evolved galaxies at z ~ 1 .5 — 3 (Daddi 
et al. 2005a; Trujillo et al. 2006; Zirm et al. 2007; Tacconi 
et al. 2006; 2008; Cimatti et al. 2008) and have high star 
formation efficiencies compared to typical galaxies of simi- 
lar mass (M - 10*" - lO'^M©; Daddi et al. 2008; see also 
Bouche et al. 2007). Their space density evolution is a cru- 
cial issue for understanding the formation of massive galaxies, 
as the SMGs are thought to represent massive galaxy mergers 
which will rapidly evolve into passive systems (Daddi et al. 
2007ab; Tacconi et al. 2008; Cimatti et al. 2008). Therefore, 
searching for the very high redshift tail of SMGs holds clues 
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into understanding early type galaxy formation and character- 
izing the distribution of their formation redshifts, both crucial 
issues for galaxy formation models in a ACDM universe. 

The strong decline of SMG counts at redshifts higher than 
z ^ 2 might signal that z ^ 2 is indeed the main formation 
epoch of massive early-type galaxy systems (Daddi et al. 
2004ab), but it can be affected by selection biases, primar- 
ily the requirement of a detection in the radio, which suffers 
from positive K-correction. On the other hand, the existence 
of a population of passive systems at z ^ 1.5 — 2.5 (Spinrad 
et al. 1997; Cimatti et al. 2004; McCarthy et al. 2004; Daddi 
et al. 2005a; Saracco et al. 2005; Kriek et al. 2008) with 
large masses (0.5 — 5 x lO'^M©), a space density of order of 
10^"* Mpc"^ (^ 10 — 30% of the coiTespondent one at z = 0; 
Daddi et al. 2005a; Labbe et al. 2005; Kong et al. 2006) 
and estimated formation redshifts for their stars of z > 3 or 
even much higher (Daddi et al. 2005a; Labbe et al. 2005; 
Maraston et al. 2006; Longhetti et al. 2007; Cimatti et al. 
2008), seems to require in turn the existence of a substantial 
population of vigorous starburst galaxies at z > 3 — 4, which 
should be detectable from their strong far-infrared (FIR) emis- 
sion. For example, if this population of z ~ 1 .5 — 3 early 
type galaxies were mainly formed during 100 Myr long bursts 
within 3 < z < 6, one would expect to find 10^^ Mpc^^ 
galaxies with SFR ^ lOOOMg, yr^^ in the same redshift range, 
which would be detectable as SMGs and with a sky density of 
« 200 deg-2. 

Evidence has been mounting in the past years that such ob- 
jects actually exist. For example, a well known case is that 
of HDF850.1 (Dunlop et al. 2004), an SMG in the Hubble 
Deep Field that is thought to be at very high redshift, ly- 
ing behind a bright foreground galaxy. In addition, although 
largely lacking spectroscopic identifications, the survey of 
Dannerbauer et al. (2002; 2004) presented evidence that a 
high fraction of sources selected with the 1 .2mm bolometer 
MAMBO at the IRAM 30 m telescope (Kreysa et al. 1998) 
are at high redshifts, given their very faint /T-band and radio 
detections or upper limits. Similarly, Younger et al. (2007) 
recently reported on seven sources selected with the AzTEC 
camera at 1.1mm (Wilson et al. 2008) with flux densities 
above 5 mJy in the COSMOS field (Scovifle et al. 2007), 
of which five (70%) are found to be very faint at optical and 
radio wavelengths, suggesting redshifts z 3 — 4. More re- 
cently, Wang et al. (2007; 2008) and Dannerbauer et al. 
(2008) presented high resolution submm/mm imaging of a 
bright SMG from GOODS-N cafled GN10'°, suggesting that 
this source is at a much higher redshift than thought before, 
probably 3.5 < z < 6. Follow-up surveys of high redshift ra- 
dio galaxies also resulted in submm detections (Archibald et 
al. 2001; Greve et al. 2007). Knudsen et al. (2006) proposed 
az = 4.048 identification for SMM J 16359+66 130. However, 
this source lacks both radio and CO identifications, and to our 
knowledge the result was never published in a refereed jour- 
nal. Recently, Capak et al. (2008) also suggest a z = 4.547 
identification for an SMG discovered with AzTEC and IRAM 
in the COSMOS field. This high redshift galaxy is > 1" away 
from the radio and IRAC positions. Offsets of up to ~ 1" are 
indeed sometimes found for SMG counterparts (Chapman et 
al. 2005). In a few of such cases the identification turned out 
to be wrong (see, e.g.. Pope et al. 2008), however Schinnerer 
et al. (2008) confirmed the z = 4.5 redshift for this galaxy. 

Here and in the remaining of the paper we adopt the "GN" nomenclature 
for GOODS-N SMGs taken from Pope et al. (2006). 



Most of the evidence for z > 4 SMGs comes from survey 
in the millimeter regime. At first glance, it is not too clear 
how to reconcile these results with much of the work done 
on SCUBA surveys follow-up. The slight difference in se- 
lection wavelength (850/im versus 1.1mm) is likely not suf- 
ficient to justify the diverging results. A possible explana- 
tion is the trend for brighter SMGs to be systematically at 
higher redshifts, as claimed by Ivison et al. (2002) and Pope 
et al. (2006; P06 hereinafter). Difficulties in locating the 
counterparts within the large submm/mm positional uncer- 
tainties, and the fact that the most distant sources could be 
ultra-faint or completely undetected in the optical, near-IR, 
mid-IR, and/or radio, imply that it is likely to obtain incorrect 
identifications by choosing the brightest galaxy within the er- 
ror circle of the submm position. Given that star formation 
proceeds through the consumption of the molecular gas reser- 
voirs, the detection of CO emission is still the most secure 
means to measure the redshift of an SMG. 

If SMGs are the progenitors of very massive early type 
galaxies, signaling perhaps their last and strongest starburst 
before the onset of passive evolution, they should also be 
strongly clustered and reside in overdense regions of the uni- 
verse. Blain et al. (2004) report evidence for significant ex- 
cess of redshift pairs in SMG samples, suggesting significant 
clustering in this population. Stevens et al. (2003) find sig- 
nificant excess of SMGs in the close surroundings of high 
redshift radio galaxy fields, evidence that SMGs might trace 
proto-cluster environments and be related to the formation of 
present-day cluster ellipticals. However, this result is not yet 
spectroscopically confirmed. More recently. Chapman et al. 
(2008a) discuss a z = 1.99 concentration of five SMGs in the 
GOODS-N field that does not correspond to a similarly strong 
overdensity of UV selected galaxies. They suggest that con- 
centration of SMGs might correspond to the cosmic structures 
with the strongest enhancement in the frequency of major 
mergers. 

In this paper, we take advantage of Plateau de Bure Inter- 
ferometer (PdBI) observations at millimeter wavelengths in 
the GOODS-N area that serendipitously revealed two strong 
emission lines, coincident with the positions of GN20, the 
brightest SMG in the GOODS-N field (^gso^m = 20.3mJy), 
and GN20.2 (^gso^^m = 9.9mJy). Using multiwavelength pho- 
tometric imaging data and Keck spectroscopy, we identify 
the lines as CO[4-3] transitions at z = 4.055 ± 0.001 and 
z = 4.051 ±0.003. Thanks to the serendipitous detection of 
CO lines in the PdBI field of view, this reverses the typi- 
cal SMG counterpart identification process in which a plausi- 
ble optical/near-IR counterpart with a measured spectroscopic 
redshift is only later confirmed through CO observations. In 
addition, we provide spectroscopic (CO) and photometric ev- 
idence for a second plausible counterpart to GN20.2, likely at 
the same redshift. This concentration of submm galaxies ap- 
pears to define a proto-cluster structure as traced by B-band 
dropout Lyman break galaxies and massive IRAC selected 
z > 3.5 galaxies. The newly measured z > 4 CO redshifts 
for SMGs have been exploited to calibrate and test the use of 
radio-IR photometric redshifts for the SMG population, and 
to obtain a new estimate of the abundance of very high red- 
shift SMGs. 

The paper is organized as follows. In Sect. |2] we describe 
the PdBI CO line detections, and determine their redshifts in 
Sect. [3] The physical properties of the sources, in terms of 
molecular gas, SFRs and masses, are discussed in Sect. |4] In 
Sect|5]we show that these objects lie at the center of a proto- 
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Fig. 1. — Velocity averaged spatial maps at 91.4 GHz centered on the serendipitously detected emission from GN20 (the 2" wide cross refers to the VLA 
1.4 GHz position). The two top panels show B-configuration 1.3" beam maps (20" field of view) while the bottom panels show the D-configuration 5.5" beam 
ones (40" field of view). Red colors correspond to positive fluxes and gray colors to negative fluxes. Contours level are shown in steps of Icr with positive 
(negative) contours shown as solid (dashed) lines. The CO[4-3] maps are averaged over 750 km s"' corresponding to the observed range of CO emission (see 
Fig. [3) and were cleaned. The continuum maps show the average of all channels outside the line range (for a 2250 km/s total bandpass). Notice that GN20 is 
observed off from the phase center with primary beam attenuation (PBA, labeled in each image) of 1.33 and 1.77 at its position in the B- and D-configuration 
observations (these maps were not corrected for attenuation in order to preserve the flatness of the noise). Given the relatively low S/N ratios, the continuum 
maps were not cleaned as the side lobes are below the noise level. 



cluster environment. We describe a method to obtain accu- 
rate photometric redshifts for SMGs in Sect|6l and discuss the 
properties and present a new selection technique for identify- 
ing high redshift starbursts in Sect]?] Throughout the paper 
magnitudes are expressed in the AB scale unless stated oth- 
erwise, and a standard cosmology is adopted with ilA,i^M = 
0.76,0.24, and h = //o/(100 km s^^ Mpc-')= 0.73. 

2. PDBI OBSERVATIONS 

We have used the IRAM PdBI to map the CO[2- 
1] transition from the z = 1.522 galaxy BzK-21000 
(J123710.60H-622234.6) redshifted to about 91.4 GHz. We 
obtained 7.5 hours of on-source integration in D-configuration 
(5.5" synthesized beam) during April 2007 and 6.4 hours on 
source follow-up observations were obtained in January 2008 



with a B-configuration (1.3" synthesized beam). For both D- 
and B-configuration observations the PdBI was used with all 
6 antennas available (15 independent baselines). The correla- 
tor consists of 8 independent units each of which was cover- 
ing 320 MHz (128 channels each with a width of 2.5 MHz) 
with a single polarization, covering a total bandwidth of about 
1 GHz with both polarizations. Observations were tuned at a 
central frequency of 91.375 GHz (for a total velocity band- 
pass of about 3000 km/s). For calibration of the data we ob- 
served standard bandpass calibrators (J0418+380, 3C273), 
phase/amplitude calibrators (J1044+719, Jl 150+497), and 
flux calibrators (MWC 349, 3C273) to which we regularly 
switched during the primary target observations. We re- 
duced the data with the GILDAS software packages CLIC 
and MAP. The maps of the fields obtained using natural 
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Fig. 2. — The same of Fig.[T]but for the galaxy GN20.2a (the 2" wide cross refers to the VLA 1.4 GHz position). The CO[4-3] maps are averaged here over 
1125 km s^' coiTesponding to the observed range of CO emission for GN20.2a (see Fig. [5). The continuum maps show the average of all channels outside the 
line range (for a 1875 km/s total bandpass). GN20.2a is even more offset than GN20 from the phase center with primary beam attenuation (PBA, labeled in each 
image) of 1.71 and 2.91 at its position in the B- and D-configuration observations. Before creating these maps, the bright emission at the position of GN20 was 
fitted in the uv data with a point source model and subtracted, in order to avoid contamination at the position of GN20.2 from the side lobes of the more luminous 
GN20 galaxy. These maps were not cleaned as the side lobes are below the noise level in all cases. No correction for PBA is applied to these maps, in order to 
preserve the flatness of the noise. 



weights and using the full 1 GHz bandpass have noise lev- 
els of about 67.9/iJy/beam and 51.2/iJy/beam for the D- and 
B-configuration data, respectively. 

Daddi et al. (2008) presented the CO[2-l] detection of the 
BzK galaxy with the D-configuration data. In addition to the 
BzK galaxy, from the initial D-configuration observations we 
also found a strong, serendipitous source in the field, located 
22" South and 8" East of the primary target (Fig. [U. Even 
when averaging all available channels over the 3000 km s^' 
spectral range observed, and after subtracting from the data 
the signal from the primary target, the serendipitous source 
was securely detected, despite the large primary beam attenu- 
ation at its distance from phase center (the primary beam re- 
sults in an approximately Gaussian response with full width at 



half maximum [FWHM] of 55". 2 at this frequency). For this 
reason, for the follow-up observations in the more extended 
B-configuration we offset the pointing by 10" from the BzK 
galaxy in order to reduce the primary beam attenuation at the 
position of the serendipitous detection without compromising 
the data quality at the position of the primary target. In addi- 
tion to confirming this serendipitous source, the new observa- 
tions also revealed a second source at high significance (Fig. [2] 
and[3]l, and a third tentative detection. The two secure detec- 
tions correspond well to the positions of GN20 and GN20.2 
in the sample of SMGs published in P06. In both cases, it is 
clear that most of the detected signal comes from an emission 
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Fig. 3. — One dimensional spectra of CO[4-3] from GN20 (left, sampled at 25 km s^') and GN20.2a (right, sampled at 75 km s^'). The flux densities have 
been corrected for primary beam attenuation. The red curve shows the best fitting Gaussian to the measured signal, allowing for the presence of faint underlying 
continuum. 



Table 1 . Radio and CO coordinates of the sources 



ID 


RAco(J2000) 


DECco(J2000) 


RAvla(J2000) 


DECvla(J2000) 


GN20 


12:37:11.90 


62:22:12.1 


12:37:11.89 


62:22:11.8 


GN20.2a 


12:37:08.77 


62:22:01.7 


12:37:08.78 


62:22:01.8 


GN20.2b 






12:37:09.73 


62:22:02.6 



Note. — The CO coordinates are from our 1.3" beam, B-configuration observations of CO[4-3]. The VLA coordinates are from the 1.4 GHz radio continuum map of Morrison et 
al. (in preparation) with 1.7" beam. 



line, most likely from CO molecular gas . Measured quanti- 
ties are summarized in Tables [T] and |2] and derived quantities 
(like luminosities and masses) are in Table [3j 

2.1. CO detection of GN20 

GN20 is very securely detected in both D- and B- 
configuration observations (Fig. [T| left). IVlost of the signal 
appears to be due to a CO emission line. Averaging the signal 
from 150 km s^' to 900 km s^' (here and in the following 
all velocities are defined with respect to the central frequency 
of the observations at 91.375 GHz) results in a detection with 
S/N = 10 for the D-configuration data (Fig.[T]left and Fig. [3] 
top). The B configuration observations within the same veloc- 
ity range results in a detection with S/N —13 and a position of 
a(y2000) = 12''37'"ir .90, 5(72000) = 62''22'"12M (Fig.[I] 
left). This position is consistent with the 7" submm error box 
of the SMG GN20 (P06), and is within 0.5" of the accurate in- 
terferometric lono et al. (2006) Sub-Millimeter- Array (SMA) 
detection of this source and of the radio counterpart in 20cm 
VLA data (Morrison et al., in preparation; Fig. [Hi. The PdBI 
position coincides well with a faint galaxy that P06 already 
identified with the optical counterpart of GN20 (Fig|5]l. GN20 
is the brightest SMG in the GOODS-N area. The radio, CO, 
and SMA submm positions actually appear to be slightly off- 
set to the East of the source seen in the optical (by ^ 0.5" or 

" It is unlikely that the emission line is not from CO. given that CO lines 
are the brightest lines observed at millimeter wavelengths. 



so), as also noted by lono et al. (2006). It is not clear if this 
is due to some additional galactic structure being present but 
extremely obscured, or just a random fluctuation across the 
data sets. 

The PdBI D- and B-configuration datasets were indepen- 
dently fitted with a point source and the resulting spectra were 
corrected for primary beam attenuation (PBA) and coadded 
with appropriate weighting (the D-configuration observations 
are noisier on an absolute flux scale, also due to the larger cor- 
rection for PBA). We performed Gaussian fitting of the result- 
ing spectrum (Fig.|3]l to determine the CO line properties, al- 
lowing for the presence of a faint underlying continuum. The 
best fitting Gaussian line is offset by 580±40 km s^' rela- 
tive to the central frequency, and has an velocity integrated 
flux of Ico = 1-5 ±0.2 Jy km s^\ with both errors given at 
the 1(7 level and following Avni (1976) for the case of two 
parameters of interest (line velocity and integrated emission 
line flux). There is evidence for faint 3mm continuum as 
well, as shown in Fig. [T] (right panels) and Fig. |3] To con- 
firm this result, we averaged the signal in UV space for veloc- 
ities outside a 1200 km s^' range centered on the line (about 
1.7xFWHM), where negligible line contribution should be 
present. The resulting continuum map has a faint source at a 
position consistent within the errors with that of the emission 
line and corresponding to a signal of 0.33 mJy, significant at 
the 5ct level. To our knowledge, this is the first secure dust 
continuum detection of an SMG in the 3mm band. 
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Fig. 4. — VLA 1.4 GHz radio continuum map of the GN20 (left) and GN20.2 (riglit) region. Tlie VLA observations were mostly taken in A-array and have a 
beam of 1.7" FWHM (Morrison et al., in preparation). The 7" radii circles are centered on the SCUBA positions from P06. The radio sources within the S50fj,m 
beam are labeled. The faint radio emission about 1" North of the GN20.2 SCUBA beam is real and corresponds to a galaxy with Zpiun = 1-6. From its mid-lR 
emission we infer only moderate SFR levels that should not be substantially contributing to the flux density observed at submm wavelength at the position of 
GN20.2. 



2.2. CO detection of GN20.2a 

Averaging the velocity channels between -300 km s^' and 
825kms-' a source is detected in the PdBI map with S /N — 
5.8 at the position of a(72000) = 12/'37'"08^77, 5(72000) = 
62''22"'01■^7 in the B -configuration observations (Fig. |2] left 
and Fig. [3] right). Positive flux with S/N just below 3 is ob- 
served also in the D-configuration data at this position (Fig.|2] 
left). We extracted the spectra by fitting point sources in the 
D- and B-configuration datasets independently and coadded 
the spectra with weighting after correcting for primary beam 
attenuation. Averaging the total resulting spectrum within - 
300 km s^' and 825 km s^^ we find a positive signal signifi- 
cant at the 7(7 level. 

The position of this source is within 0.2" of the relatively 
bright radio source that P06 identified with the counterpart of 
the SMG GN20.2, close but well distinguished from GN20 
that is 24" to the NE. The position corresponds to a faint 
galaxy in the deep Hubble Space Telescope {HST) images 
(Fig. |5]). Also in this case the signal appears to be mostly 
due to an emission line. The best fitting Gaussian line is off- 
set by 330±150 km s^' relative to the central frequency, and 
has a velocity integrated flux of Ico = 0.9±0.3 Jy km s^'. 
The largest contribution to the flux error is the uncertainty in 
the continuum level. Some positive continuum appears to be 
present, although the S/N is not yet sufficient for a reliable 
measurement (Tab. |2]i. Accounting for the possible contin- 
uum emission and its uncertainty, the CO emission line is still 
detected at the 6a level. 

2.3. GN20.2b: a second counterpart to GN20.2 ? 

Inspection of the newly reprocessed 1.4 GHz radio map 
around the position of GN20.2 (Fig. 21 Morrison et al., in 
preparation) shows an additional faint radio source with a 
flux density of 32.2^Jy (5cr) at a(72000) = 12''37"'09^73, 
(5(72000) = 62''22'"02".6 (radio frame), 3.4" from the 



GN20.2 SCUBA position (P06) and 6.8" from the CO de- 
tected counterpart discussed above. We refer to the primary 
counterpart of GN20.2 as GN20.2a, and this newly identified 
faint radio source is referred to as GN20.2b. GN20.2b is a 
plausible additional counterpart to the submm emission. 

Interestingly, positive signal is found in the PdBI map close 
to the position of GN20.2b. The PdBI D- and B-configuration 
datasets were independently fit with a point source fixed at 
the VLA position of GN20.2b (a secure position cannot be 
derived from the PdBI data alone, due to the low S/N) and 
the resulting spectra were corrected for primary beam atten- 
uation (PBA) and coadded weighting according to the noise, 
similarly to what was done for the other galaxies. Averaging 
over the whole 1 GHz bandwidth in the coadded spectrum re- 
sults in a ^ 3(7 signal. If we average the spectrum between 
« km s^' and w 700 km s^', we find a > 4(7 signal. Al- 
though longer integrations at PdBI would be required for a 
more reliable conclusion, the data are again consistent with 
CO emission, quite close in frequency to the other securely 
detected lines. If real, that would correspond to an integrated 
flux of Ico ~ 0.45 Jy km s"'. 

3. REDSHIFT IDENTIFICATIONS 

Depending on which transition we are actually observing, 
the CO detection fixes the redshift of GN20 and GN20.2a 
to 1 + zco ^ 1 -264 X n, where the observed transition is 
CO[{n) - (n - 1)]. Possible redshifts ai-e: 0.26, 1.53, 2.80, 
4.05, 5.32, 6.58, etc. The most straightforward guess for the 
redshift would be 1.53, corresponding to companion galax- 
ies to BzK-21000, the primary target of the PdBI observa- 
tions. However, we show in the following that the correct 
identification for the lines is CO[4-3] at redshift zco ^ 4.05, 
making these among the most distant CO detections so far 
for SMGs. Although it might seem unlikely that we have 
serendipitously detected CO emission lines from z ^ 4 galax- 
ies, we note that the 1 GHz spectral range of our data maps 
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Fig. 5. — Multiband HST+ACS imaging of the optical counterparts of GN20 (top), GN20.2a (second from top), GN20.2b (third from top) and a B-dropout 
Lyman break galaxy in the area having also Zspe^ = 4.058 from Keck spectroscopy (the coordinates of this source, named BD29079 in the GOODS-N B-dropout 
catalog, are 0(72000) = 12''37'" 1 1'.53, <5(/2000) = 62''21"'55'.7). From left to right, as labeled, the data refer to the ACS filters F435W, F606W, F775W and 
F850LP (sometimes referred to in the text also as B-, v-, i- and z-band, respectively). The images have been smoothed with a 3 pixel FWHM Gaussian. Circles 
around the sources have 1" radius and are simply centered on the identified ACS counterparts. 



into a twice larger redshift range at z 4 versus z ^ 1 .5, which 
in turn corresponds to a two times larger comoving volume. 
Also, at higher redshift the density of CO transitions per unit 
frequency and at fixed observed frequency is higher, hence 
chance detections are non-negligible once a luminous high 
redshift target is within the field of view (Fig.|6]l. 

The following section discusses several independent red- 
shift estimates for these galaxies, e.g., based on optical pho- 
tometric redshifts, radio-IR photometric redshifts, and optical 
spectroscopy. This comparison is interesting as the problem 
of unidentified CO line detections will become a typical situa- 
tion of the coming years with new, wide-band instrumentation 
rapidly becoming available. Even 8 — 20 GHz receivers will 
often identify single CO lines from single galaxies (Fig. |6]l. 
A direct assessment and inter comparison of multiple redshift 
identification techniques is thus interesting for this first test 



case of serendipitously detected CO lines. 

3.1. GN20 

3.1.1. Photometric redshift from stellar emission 

Fig. H] shows HST+ ACS imaging in the four bands avail- 
able in GOODS (Giavalisco et al. 2004a). As already noted 
by P06, the optical counterpart of GN20 is a B-dropout Ly- 
man break galaxy according to the definition of Giavalisco 
et al. (2004b). This definition returns galaxy samples with 
z = 3.78 ±0.34, favoring the identification of the PdBI line 
as CO[4-3] at z = 4.055. Circumstantial evidence for the 
likelihood of a redshift around z = 4.05 had been even ear- 
lier presumed from the presence of a bright B-band dropout 
galaxy at a 16" separation (Fig.|5]l with fairly similar colors 
and a known spectroscopic redshift of z = 4.058 measured 
from Keck spectroscopy (Stern et al., in preparation, we refer 
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Redshift 

Fig. 6. — The frequency-redshift relation for CO emission lines in the 
3mm atmospheric window. The dotted range around 91.4 GHz corresponds 
to the PdBl observations presented in this paper. The dashed line centered at 
100 GHz coiTesponds to an 8 GHz range, as planned for future generations 
of correlators, like those that will be used on ALMA. 



to this object as BD29079). The IRAC colors of the GN20 
Z?-band dropout counterpart are very red, with the brightest 
IRAC band being the S.O^m (see, e.g., P06). If due to the 
emission of stars, this implies that the 1.6/imrest frame bump 
is beyond the 5.8/im channel, or z 3.5. However, the IRAC 
emission can also be red due to, or affected by, the presence of 
obscured AGNs (e.g.. Stern et al. 2005; Daddi et al. 2007b). 

Stellar photometric redshifts were thus estimated using only 
the four HST bands and running the hyperz code (Bolzonella 
et al. 2002). The Maraston (2005) models were used, with a 
variety of star formation histories and allowing for reddening 
using the Calzetti et al. (2000) law. The near-IR bands are not 
used because the available datasets, described, e.g., in Daddi 
et al. (2007a), are not sufficiently deep for placing meaning- 
ful constraints on the spectral energy distribution SED of this 
faint and relatively blue galaxy. The distribution of values 
versus redshift (Fig|7]i is consistent with 3.4 < z < 4.2 (95% 
confidence range). Fig. [T] shows the observed SED and the 
best fitting model for the case of z = 4.055. 

3. 1 .2. Photometric redshift from the infrared and radio 
emission 

An independent estimate of the redshift of GN20 can be de- 
rived from its FIR emission. The ratio between mid-IR, FIR, 
mm and radio flux densities is redshift dependent. This is 
due to the fact that the radio luminosity is observed locally 
to be proportional to the bolometric IR luminosity (Condon 
et al. 1992; Yun et al. 2001), which seems to hold at high 
redshift as well (Elbaz et al. 2002; GaiTett 2002; Appleton 
et al. 2004). Originally this idea was exploited by Carilli 
& Yun (1999; 2000) to estimate photometric redshifts of dis- 
tant SMGs, also attempted by several other studies (Hughes 
et al. 2002; Wiklind 2003; Aretxaga et al. 2003; 2005; 2007; 
Clements et al. 2008). In addition, one can exploit the fact that 



the mid to FIR emission has a generally well defined peak at 
50 — 200/im, affecting directly the ratio of mid-IR to submm 
and mm emission. 

For this purpose, we use the available photometry at 24/Ltm, 
850/im, 3mm and 20cm. Some of the Spitzer+MWS 24^m 
and VLA 20cm measurements are updated from P06 using the 
most up-to-date GOODS datasets. We measure a 24-fim flux 
density of 65.5 ±3.5 /iJy and a 20cm flux density of 75.8 ± 
7.9 /iJy. We use a 850/im flux density of 20.3 ±2.3 mJy from 
P06 and a 3.3mm flux density of 0.33 ±0.07 mJy (Section|2li. 

The observed flux densities have been compared, as a func- 
tion of redshift, to the predictions of a suite of 105 tem- 
plate SEDs that were built following the luminosity corre- 
lations observed for local galaxies as described in Chary & 
Elbaz (2001). The template dust temperatures get warmer 
with increasing luminosity from a minimum template lumi- 
nosity of Log Lir/L© = 8.41 up to a maximum luminosity 
of Log Lir/Lq = 13.55. Radio continuum emission is added 
to each template following the radio-IR correlation (Yun et 
al. 2001). For each template of a given luminosity we com- 
pute the expected flux densities in the available bands (24/im, 
850/im, 3mm and 20cm) as a function of redshift and per- 
form a minimization over the template total IR luminos- 
ity as a function of redshift, without allowing template nor- 
malizations to vary. In order to avoid biasing the fitting to 
a particular band, we increased the formal uncertainty in the 
24/im flux density to obtain a S/N ratio of 9 (as opposed to 
the higher measured value of ^ 20), comparable to the radio 
and submm bands. 

FiglS] shows the results. The best fitting redshift is Zphot = 
3.7 and at the 99% confidence level is 3.3 < z < 4.3. There- 
fore, we conclude that the radio-IR SED also independently 
supports the z — 4.055 redshift determination. We notice that 
this analysis allows us to reject with high confidence any red- 
shift identification with z < 3 or z > 5. 

FigS (right panel) also shows the comparison of the ob- 
served GN20 SED from UV to radio to the best fitting (and 
brightest) template in the library of Chary & Elbaz (2001). 
The SED is remarkably well reproduced in the mid- to FIR 
and radio, without any rescaling or ad-hoc normalizations. 
The UV-optical part of the template SED also happens to 
match the shape of the observed data quite well, although 
there is a factor of ~ 2 offset over the IRAC bands that im- 
plies a different L(FIR)/L(optical) ratio, suggesting a higher 
specific SFR (SSFR). A similarly good agreement, including 
the optical/UV part, was seen by Daddi et al. (2005b) when 
comparing the average SED of z = 2 BzK selected ULIRGs 
to Chary & Elbaz (2001) templates. The optical/UV emis- 
sion was simply added to the models using a synthetic stellar 
population model matched to the Arp220 SED and scaled to 
reproduce the local broad correlation between bolometric lu- 
minosity and B-band rest frame luminosity. 

3.1.3. Keck spectroscopy 

A solid identification of the CO line, and thus a deter- 
mination of GN20's redshift, is finally provided by deep 
optical spectroscopy. The counterpart of GN20 was ob- 
served with the Deep Imaging Multi-Object Spectrograph 
(DEIMOS; Faber et al. 2003) on the Keck II telescope several 
times between 2004 and 2007, though only the final observa- 
tions on UT 2007 April 14—15 were in good conditions. The 
total integration time was 2.5 hr, and a 1"2 wide slit was ori- 
ented along the major axis of the galaxy (PA = 13.8°). The 
600ZD grating (A^iaze — 7500 A) was used, providing a res- 
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al. (2006) and we show an upper limit of 3 mjy. 



olution R = 2000. The GG455 order-blocking filter was in- 
stalled. We processed the data using a slightly modified ver- 
sion of the DEEP2 DEIMOS pipeline^^ 

The data show no emission lines while the continuum has 
a low S/N per pixel, as expected given the faintness of the 
source (/ ~ 25) and its spatial extent of about 1.5", implying 
a low surface brightness. The emission is nevertheless clearly 
seen after smoothing the spectrum. We extracted the galaxy 



See |http : //astron . berkeley ■ edu/~cooper/deep/spec2d/ 



spectrum with apall under IRAF using a 1.2" wide aperture. 
We measured the rm.s. per pixel as a function of wavelength 
based on the sky spectrum. This information was used to 
properly weight pixels when binning by wavelength (e.g., to 
down weight noisy pixels affected by OH sky lines) and to ac- 
curately estimate the error of the binned spectrum. Given that 
for z = 4.055 the Lya forest break would fall at 6 145 A we in- 
cluded two spectral bins bracketing this wavelength to verify 
if a break is indeed present as expected. Fig. [T] (right panel) 
shows the results. Although the Keck spectrum only extends 
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Table 2. Measured properties of GN20, GN20.2A and GN20.2b 



ID 


^co 


ZKeck 




5'3.3mm 


Sl.4GHz 








4:014-31 


FWHM 








(mJy) 


(mJy) 


(MJy) 


(MJy) 


(MJy) 


(MJy) 


(Jy km s^') 


(kms-') 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


GN20 


4.055 ±0.001 


4.06±0.02 


20.3±2.lt 


0.33±0.07 


75.8±7.9 


65.5±3.5 


26.1 ±2.6 


9.5±1.0 


1.5±0.2 


710±120 


GN20.2a 


4.051 ±0.003 


4.059 ±0.007 


< 9.9±2.3t 


<0.20 


180.7±8.4+ 


30.2±5.6 


9.8±1.0 


4.1 ±0.4 


0.9 ±0.3 


1100 ±400 


GN20.2b 






<9.9±2.3+ 


<0.18 


32.2±6.5 


12.0 ±4.3 


15.4±1.5 


10.1 ±1.0 







Note. — Col. (1): Name of the object. Col. (2): Redshift derived from the CO[4-3] observations. Col. (3): Redshift derived from the Keclt optical spectroscopy. Col. (4): SCUBA 
850 fim flux density from P06^. Col. (5): Continuum flux density at 3.3mm. We notice that only GN20 is detected in the 3.3mm continuum. For GN20.2a and GN20.2b we list 2-sigma 
upper limits. The formal 3.3mm continuum measurement are: 0. 1 1 ± 0. 10 and 0.08 ± 0.09, respectively. Col. (6): VLA 1.4 GHz flux density from P06^ and our own measurements. 
Col. (7): Spilzer+MlPS 24 ^im flux density. Col. (8) and (9): Spitzer+IRAC flux densities. All detections have high significance, eiTors are dominated by systematics effects at the 10% 
level. Col. (10): Velocity integrated CO[4-3] flux. Col. (11): Width of CO[4-3] line, estimated through a single Gaussian fitting. 



down to 5900A, and therefore only one binning point is avail- 
able with decent S/N below the expected Lya forest break, 
we find that the spectrum strongly supports the z — 4.055 red- 
shift and in particular it reveals the presence of the expected 
Lya forest break. The flux density measurement blueward of 
the expected break corresponds to a Lya forest decrement of 
Da = 0.55 ±0.10 (derived using the stacked bins closer to the 
break), fully consistent to what is expected for the mean trans- 
mission of the intergalactic medium at z = 4.05 (Songaila et 
al. 2004). When running a photometric redshift analysis on 
the binned Keck+DEIMOS data (Fig.|2] left panel) we see that 
of the possible CO line transition identifications, the z = 4.055 
redshift is unambiguously chosen with very high confidence. 

3.2. GN20.2a 

A nearby (24" separation; 169 kpc proper at z = 4.05) SMG 
companion to GN20 is GN20.2, first discovered by Chap- 
man et al. (2001a)'l GN20.2 has an 850^m flux density 
of 9.9 mJy as reported by P06, that also identifies the most 
likely optical counterpart with a relatively bright radio source 
(5'(20cm) = 180 fiiy). This galaxy is also a B-drop (P06; 
Fig.|5]l with a faint 24/im flux density of 30.2/iJy. 

An analysis similar to the one performed in Sect l3.1 [ results 
in a photometric redshift of Zphot.opt = 3.93 and a 90% con- 
fidence range within 3.55 < z < 4.05. Similarly, the analy- 
sis of the IR SED suggests a redshift of 3.4 < z < 4.2 at the 
99% confidence. This confirms that the overall properties of 
this source are compatible with being at the same redshift of 
GN20. This is further demonstrated by deep Keck+DEIMOS 
spectroscopy (2.5 hours) of this galaxy that was obtained in 
April 2005 during cirrus conditions, and reduc ed in the same 
way as for the GN20 observations (see Sect. 13.1.3b . Fig. |9] 
shows that a strong spectral break is detected at about 61 50A, 
consistent with the location expected for the Lya forest break 
at a redshift of 4.055. The observed break corresponds to 
a Lya forest decrement Da ~ 0.67 ± 0.05, larger than the 
average expectations but still within the observed scatter in 
the intergalactic medium transmission at z = 4.05 (Songaila 
et al. 2004). We constrain the redshift of this source to 
be z = 4.059 ±0.007 (la), based on fitting the Keck spec- 
troscopy with the Maraston (2005) model shown in Fig. |9] 
demonstrating that this source is at a similar redshift of its 
brighter companion GN20. We conclude that also in this case 
the detected CO emission is CO[4-3] at z = 4.051 ±0.003, in 
good agreement with the Keck spectroscopy. 

" Chapman et al. (2001a) used SCUBA to observe (and successfully de- 
tect) objects from a sample of optically faint radio sources with / >25 mag. 
GN20.2, alias VLA J123708±622201, was part of this survey. 



3.3. GN20.2b: redshift and contribution to the submmflux 

The GN20.2b radio source is 3.4" from the GN20.2 
SCUBA position of P06, closer than GN20.2a which is 6.0" 
(although, we caution, that the SCUBA position of GN20.2 
might be less accurate than for the average SMG, due to the 
brighter and nearby companion, GN20). The radio position is 
close (only 0.5" to the north) to a very faint galaxy (z = 27.34) 
detected in the HST+ACS imaging (Fig|5]l. It is possible that 
this faint radio galaxy could be at the same z ^ 4.05 redshift 
of GN20 and GN20.2a, and it could be contributing some of 
the submm flux detected at the position of GN20.2. We will 
discuss these two issues in turn in the following. 

Although undetected in the F435W band, the faint ACS 
counterpart to GN20.2b is not formally classified as a Z?-band 
dropout due to the upper limit in the F435W not being strin- 
gent enough, but its blue (/ — z) and red (v — /) colors support 
the hypothesis that also this object is at z 4. The galaxy is 
also fairly bright and red at IRAC wavelengths (SED peak at 
8.0/im with 20.93 mag) and is selected as a z > 3.5 massive 
galaxy candidate in the work of Mancini et al. (2008). In- 
specting the highest resolution IRAC imaging at 3.6/im and 
4.5/im it appears that blending with a neighbouring source is 
somewhat affecting the Spitzer IRAC measurements and bias- 
ing the flux densities high (likely by at most a factor of 1.5- 
2). We used the radio-IR photometric redshift technique to 
constrain its redshift, in the assumption that the emission is 
dominated by star formation and not an AGN. We find that by 
itself, the large radio to 24/im flux density ratio for this source 
strongly constrains its redshift to z > 3.2 (99% confidence 
level), independently of how much it contributes to the 9.9mJy 
flux density at 850/im. This is consistent with this object be- 
ing an intrinsically very luminous, highly star forming galaxy 
at high redshift (although, in principle, the high ratio could 
be also due to radio-lo ud A GN emission from lower redshift). 
As discussed in Sect. 12.31 there is possible evidence for CO 
emission from this source, where a positive signal is found at 
the > 4a level. Although the situation is not secure as for the 
other two sources, we conclude that there is some support- 
ing evidence for a similar z ^ 4.05 redshift also for GN20.2b. 
Taking the CO signal at face value, we would assign a redshift 
of z = 4.052 ±0.006 to GN20.2b, where the error corresponds 
to half the frequency range spanned by the tentative posi- 
tive emission, if we identify that as CO[4-3]. This is within 
Az = 0.003 of the redshift of GN20 and GN20.2a, while our 
data is sensitive in principle to emission from a Az = 0.05, 17 
times larger redshift range of 4.02 < z < 4.07. We conclude 
that it is at least plausible that GN20.2b could be lying at a 
similar redshift of GN20.2a. In any case, this appear to be a 
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Fig. 9. — Keck spectroscopy results for the GN20.2a optical counterpart. The left panel shows the Keck spectrum, smoothed with a 17A boxcar (27 pixels) in 
order to increase the S/N. The Lya forest break is well detected, yielding a redshift of z = 4.059 ± 0.007 (Ict, determined from spectral fitting with a Maraston 
[2005; M05] template). This is consistent with z = 4.051 ±0.003 derived from the CO emission. The relative sky noise versus wavelength is shown at the bottom 
of the plot. The right panel is as in Fig. |7] but for the GN20.2a galaxy, showing the binned Keck spectra and the HST photometry. The best fitting M05 model is 
shown, redshifted to ^ = 4.059. 



fairly high redshift {z > 3.5) galaxy. 

Could GN20.2b be contributing some major fraction of 
the submm flux of GN20.2 ? Finding multi-component sys- 
tems is not a new situation for SMG identification — e.g., 
SMMJ094303 (Tacconi et al. 2006) and GN19 (HDF242; 
Tacconi et al. 2008) are two cases where two-component 
systems, both radio detected, have been confirmed through 
CO observations. Dannerbauer et al. (2004), P06, Ivison et 
al. (2007), and Younger et al. (2008a) report other cases of 
submm/mm galaxies with two radio counterparts. If this is 
also happening in the present case of GN20.2 is not easy to 
assess quantitatively. We present below circumstantial evi- 
dences in favor or against this hypothesis. 

The main supporting evidence is that its VLA 1 .4 GHz radio 
continuum flux density is 42 ± 10% that of GN20, which is 
similar to the ratio between the submm fluxes of the same 
two galaxies, suggesting that GN20.2b might be contributing 
some substantial flux, unless the radio emission is affected by 
an AGN. We notice though that GN20.2b is 5.5 times fainter 
in the VLA 1.4 GHz radio continuum than GN20.2a, but this 
is likely due to the fact that the radio emission of GN20.2a 
(even more than twice larger than that of GN20) is powered 
by an AGN, as discussed in the following section. 

Evidences disfavoring instead a major submm flux contri- 
bution are the following. Comparing the 24/im flux density 
measurements, we find that GN20.2b (only tentatively de- 
tected at 24/im) is at least 2-3 times fainter at 24/im than 
GN20.2a and at least 5.5 times fainter than GN20. Account- 
ing also for blending, the 5(24/im)/5(4.5/im) flux density 
ratios are at most a factor of 2 for GN20.2b and about fac- 
tors of 6-7 for GN20 and GN20.2a (Tab. E]). This is evi- 
dence for lower specific SFR in GN20.2b, suggesting that this 
galaxy has already passed its major peak of activity. Sim- 
ilarly, the ratio of the CO fluxes detected for GN20.2 and 
GN20.2a is roughly a factor of 2, similar to the ratio of submm 
fluxes between the two galaxies (Tab.|2]). Even if the tentative 



CO emission of GN20.2b is real, this would still be at most 
half of the CO emission of GN20.2a (but, of course, there 
could be place for stronger CO emission if the galaxy lies 
outside of 4.02 < z < 4.07). Given that the CO and bolo- 
metric luminosities of SMGs are known to correlate (Greve 
et al. 2005; Solomon & van den Bout 2005; see also Sec- 
tion m, this could suggest that GN20.2a does account for 
most of the submm emission of GN20 from P06, and that 
GN20.2b could contribute at most some 30% of it but likely 
not much more. A final argument is that GN20.2 was orig- 
inally discovered by Chapman et al (2001a) in photometry 
mode with SCUBA at JCMT, pointing at the position of the 
GN20.2a radio galaxy. This resulted in a submm flux den- 
sity of 5'(850/im) = 10.2 ±2.7 that is fully consistent with the 
flux in flie P06 map. GN20.2b is 6.8" away from GN20.2a 
and any emission from its position would have been close to 
the edge of the telescope primary beam, resulting attenuated 
by a factor of 2. Nevertheless, the Chapman et al (2001a) 
flux is actually slightly brighter than the flux in P06, although 
consistent within the errors. Accounting for the errors in the 
P06 and Chapman et al. (2001a) 850/im flux density mea- 
surements, this suggests that GN20.2a does contribute at least 
half of the flux of GN20.2, and possibly most of it. 

For completeness, we notice that an additional faint radio- 
detected galaxy is present just 1" North outside the edge of 
the SCUBA beam. This galaxy is a hard X-ray source in the 
castalog of Alexander et al. (2003), for which we derive a 
photometric redshift of Zphot = 1.6. The radio and mid-lR 
emission of this object, if entirely due to star formation would 
correspond to about 1 00 M© yr^ ' . This in turn would produce 
at most a submm flux of order of 1 mJy at 850/im. Being at 
8" from the P06 position and at 14" from the Chapman et al. 
(2001a) pointed observations, any potential submm emission 
from this galaxy would be negligible for the GN20.2 system 
850/im photometry. 

In conclusion, our results suggest a scenario in which 
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Fig. 10.— The 850^m to lAfim and 1.4 GHz flux density ratios for GN20 
and GN20.2a, filled circles. We assign to GN20.2a the full 850^m flux den- 
sity of GN20.2 although G N20 .2b might be contributing some of the submm 
flux, as discussed in Sect. 13.31 In that case, both flux density ratios should 
be taken as upper limits. The lines show the predicted ratios labeled as a 
function of bolometric luminosity (Lir) for the CEOl library of templates at 
aredshiftofz = 4.055. 



two z ^ 4.05 counterparts might be contributing to the 
submm emission of GN20.2, GN20.2a and GN20.2b, al- 
though GN20.2a is hkely accounting for most of the submm 
emission. 

4. THE PHYSICAL PROPERTIES OF z = 4.05 SMGS 

4. 1 . Bolometric IR luminosities 

Following the analysis described in Sec. 13.1.21 the radio-lR 
SED of GN20 is best reproduced with a CEOl library tem- 
plate corresponding to a luminosity of Ljr = 3.5 x W^Lq, 
redshifted to z = 4.055. However, the CEOl library implies 
a correlation between dust temperature and luminosity that is 
calibrated on observations in the local Universe but that might 
not necessarily hold at higher redshifts. In order to explore 
the uncertainties in the determination of GN20's total IR lu- 
minosity due to possible SED temperature variations, we used 
again the CEOl library, but allowing this time for a free nor- 
malization of the templates when comparing to the observed 
fluxes, thus spanning the full range of dust temperatures in 
the library. For this exercise, we used the measurements or 
upper limits at 24fim, 850/im, 3.3mm and 20cm and included 
also the 1.1mm and 1.2mm measurements from Perera et al 
(2008) and Greve et al. (2008) and the 70^m upper limit 
from Frayer et al. (2006). Integrating over the best fitting 
template yields in this way Lir = 2.4 x lO'^Lp and a formal 
uncertainty from the fit at Icr of less than 0.05 dex (the lumi- 
nosity step in the CEOl library). The best fitting template has 
an intrinsic Lir = 1 .5 x lO'^L© and is scaled up in luminosity 
by a factor of 16. This would correspond to an SED temper- 
ature of about 57 K, lower than those of CEOl templates of 
comparable luminosities Ljr >, IO'^'L©. The same effect can 
be appreciated more directly in Fig. [TO] where we show the 



850/im to 24-fim and radio flux density ratios. The GN20 flux 
density ratios are reasonably in line with CEOl predictions, 
although the 850/im to radio flux density ratio is 80% higher 
than expected for the Lir = lO'^^'^Lp! galaxy in the CEOl li- 
brary. The effect is significant at the 2-a level only, but indeed 
its SED is closer to the Lir = lO'^L© expected ratio. SMGs 
have been generally found to be somewhat colder than CEOl 
models, as a result of the selection at 850/im (P06), so this 
result is not too surprising. 

As a third estimate of the IR luminosity of GN20 we might 
use directly the radio-lR correlation (the other two estimates 
also used this correlation to fit the radio, although this was 
weighted together with all the other measurements). For its 
measured 1.4 GHz flux density of 75.8/iJy, assuming a ra- 
dio continuum {f^ °^ v^°') with a = 0.8, we would derive 
Lir = 2.8 X lO'^Lp;. It is not obvious to decide which of the 
3 estimates is more reliable but fortunately the measurements 
are very close, all within a range of less than 0.2 dex. In sum- 
mary, we adopt the average value of Lir = 2.9 x 10 '^L© for 
our best estimate and an uncertainty of 0.2 dex. 

The 850/im to VLA radio flux density ratio for GN20.2a 
is about 5x lower than that for GN20 (FigfTOl). implying a 
strong radio excess of a factor of 4 — 5. GN20.2a has a radio 
luminosity of 2 x 10^^ W Hz^' at 1.4 GHz rest frame, using 
again a radio continuum index a = 0.8. This demonstrates 
that a powerful radio-loud AGN is hidden in this z — 4.05 
source. Therefore, we cannot use the 1.4 GHz flux density 
to estimate its total IR luminosity. We used the CEOl library 
templates, with and without allowing for a free normaliza- 
tion, to fit the 24/im, 70/im, 850/im and 3.3mm measure- 
ments or upper limits. Similarly to GN20, we find Lir = 
2.2 X IO^^Lq and Lir = 1.0 x IO'^Lq, flie latter being the 
'free-normalization' case that also prefers a lower tempera- 
ture template. These estimates use the full 850/im flux den- 
sity measurements although some part of this might be due to 
GN20.2b. This also assumes that the AGN inside GN20.2a 
is not contributing substantially to the measured 24/im and 
850/im flux densities. Pope et al. 2008 showed that this 
is generally the case at least for z ~ 2-3 SMGs. Given the 
uncertainties, we adopt the average of these two estimates 
Lir = 1.6 X IO^^'L© and an uncertainty of about 0.3 dex. 

This analysis implies that GN20 and GN20.2 are extremely 
luminous galaxies, with Lir > IO'^Lq. Besides being among 
the most distant known SMGs, these are also some of the 
most luminous galaxies known so far. We use a Kennicutt 
et al. (1998) conversion of SFR\WIq yr"'] = Lir[L0]/1O"', 
expressed for the Chabrier 2003 IMF adopted in this paper 
The IR luminosities correspond to SFRs of > IOOOMq yr^' 
if the IR emission is dominated by star formation. We note 
that neither the GN20 nor GN20.2a (or GN20.2b) counter- 
parts are detected in the X-rays in the catalog of Alexander et 
al. (2003), although the radio emission does suggest the pres- 
ence of an AGN inside GN20.2a. There is no clear evidence, 
instead, for the presence of an AGN inside GN20. If a top- 
heavy IMF is adopted, the implied SFR could also become 
significantly smaller. 

4.2. CO luminosities and molecular gas masses 
The observed CO[4-3] fluxes convert to luminosities of 
^co[4-3] = 6-2 X lO'^L© and L^^j^ gj = 3.7 x 10'%, for 
GN20 and GN20.2, respectively. FigHU shows that these 
galaxies appear to lie close to the (non-linear) correlation be- 
tween Lqq and Lir traced by local ULlRGs and distant SMGs. 
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Fig. 1 1 . — Comparison of FIR to CO luminosities for a variety of local and distant galaxies (see Solomon & van den Bout 2005; Greve et al. 2005; figure 
adapted from Daddi et al. 2008) and the GN20 and GN20.2a galaxies (points with error bars). The FIR estimates for GN20 and GN20.2a are described in the 
text, based on the radio-IR SEDs. The CO luminosities are estimated from the CO[4-3] transitions for GN20 and GN20.2a, and for a variety of other transitions 
for the other sources. We assume constant brightness temperature in the various CO transitions for this comparison. Notice that CO luminosities (and therefore 
gas masses) could be underestimated if the observed CO transitions are not thermalized. 



Table 3 . Derived properties of GN20 and GN20.2a 



ID 




^CO|4-3] 


il.4 GHz 




^gas 


Mdyn 




Lr.., 


K km s"' pc^ 


WHz-' 




Mp, 


sin i'^^ Maj 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


GN20 


2.9 X 10" 


6.2 X 10'° 


8.4 X lO^'' 


2.3 X 10" 


5.0 X 10'" 


2.3 X 10" 


GN20.2a 


1.6 X 10" 


3.7 X 10'° 


2 X 10^5 


0.5 X 10" 


3.0 X 10'" 





Note. — Col. (1): Name of the object. Col. (2): total IR luminosities for GN20 and GN20.2a are derived by fitting the global IR SED using all available flux density measurements 
from 24/i to 1.4 GHz with CEOl models. For GN20.2a, the 1.4 GHz measurement was not used, as it is affected likely by an AGN, and we attributed the whole 850/im emission to 
this source. Typical eiTors are estimated to be of order of 0.2 dex and 0.3 dex for GN20 and GN20.2a, respectively. Col. (3): luminosity of the CO[4-3] transition. Col. (4): 1.4 GHz 
rest frame luminosities, estimated using a radio continuum index a — 0.8 (/^^ oc ), Lj 4ghz — 47rD^5(1.4GHz)( I -\-z)~^'~ iPi is the luminosity distance). Col. (5): stellar masses 
are estimated from SED fitting from HST+ACS to Spitzer+IRAC photometiy, using constant star formation rate models from Maraston (2005) and allowing for dust reddening with 
a Calzetti et al. (2000) law. Typical errors are found to be about 0.2 dex. Col. (6): molecular gas masses assume that the CO[1-0] and CO[4-3] transitions coixespond to the same 
brightness temperature and a conversion factor Xco = 0.8 Mj, (K km s"' pc-)"' . Col. (7): the dynamical mass for GN20 is estimated following Solomon & van den Bout (2005) using 
Mdyn X sinp- /Mq — 233.5 x v^^-^ x (with velocity expressed in km s"^ and half light radius in pc). 



In the assumption that the CO[1-0] and CO[4-3] transitions 
correspond to the same brightness temperature and assuming 
Xco ^ 0.8 Maj (K km s"' pc^)-' (Downes & Solomon 1998; 
Solomon & van den Bout 2005), the CO luminosities convert 
to total molecular gas masses of 5 X IO'^Mq and 3 x lO'^M© 
for GN20 and GN20.2a, respectively. These might be re- 
garded as lower limits as, in the case that the CO[4-3] tran- 
sitions observed were not thermalized, the total CO luminosi- 
ties and gas masses would be higher VLA observations of 
CO[1-0] will be able to address this point. 

At the observed SFR levels, the gas reservoirs would be 
exhausted in roughly 20-30 Myr (FigfTTTl. Similar timescales 
have been derived for z ^ 2-3 SMGs by Greve et al. (2005). 

4.3. Dynamical and stellar mass estimates 

The CO[4-3] emission from GN20 has a FWHM of 
710±120 km s^', matching well to the typical value of CO- 
detected SMGs at lower redshifts (Greve et al. 2005). From 
the higher resolution B configuration PdBl data alone we 
find that the CO line is marginally resolved spatially, with a 



FWHM of 0.42" ± 0.18" from Gaussian fitting of the uv visi- 
bilities performed within GILDAS. The possible extension of 
the source is supported by the fact that for point source ex- 
tractions, the flux measured with the D-configuration data is 
higher than the one from the B-configuration data, a 2a ef- 
fect suggesting that the CO emission is starting to be resolved 
at the 1.3" resolution of our B-configuration data. This size 
matches well to the typical CO sizes of r,. ~ 2 kpc for SMGs 
(Tacconi et al. 2006; 2008), and is roughly consistent with 
its UV extension as well (Fig. |5]l. Recently, Younger et al. 
(2008b) measured a size of 0.6" ± 0.15" for GN20 (Gaussian 
fitting), from the far-lR dust continuum emission at 890/im 
with the Sub Millimeter Array, which is fully consistent with 
our CO size estimate. The velocity and size correspond to a 
total dynamical mass of ^ 2.3 x lO"sin/^^M0 including gas, 
stars and dark matter. 

Given the molecular gas mass estimates, this corresponds to 
a molecular gas fraction of about 20%, similar to what is typi- 
cally found for lower redshift SMGs (Greve et al. 2005). The 
stellar mass derived from SED fitting of the ACS to IRAC 
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2 4 10^ 10'" 10" 

Redshitt Stellar Mass [M^] 

Fig. 12. — Specific .star formation rate (SSFR; SFR to stellar mass ratio) as a function of redshift is shown in the left panel. The right panel shows the stellar 
mass versus star formation rate. Large solid circles are GN20 and GN20.2a. Their stellar' masses were derived through SED fitting, and typical en'ors are found 
to be of the order of 0.2 dex. The filled square (left panel) and small crosses (right panel) refer to measurements derived for the IRAC detected, massive B-band 
dropout Lyman break galaxies. Stellar masses and SFRs of the B-band dropouts were obtained through fitting of the UV to IRAC SEDs. Individual measurements 
are shown in the right panel, while the left panel shows B-dropout's average quantities (with error bars showing the error in the mean of the SSFR of the ensemble 
and the expected range in redshift). Stacking at 850/im and 1.4 GHz yield upper limits on the average SFR of order of 55M,:: yr"' for the B-band dropouts, 
consistent with the average UV measurements. The star refers to the average of SMGs at z ~ 2-3 from Tacconi et al. (2006); error bars refer to the spread of 
properties in the sample. Filled black circles in the left panel show the average SSFR versus redshift at a stellar mass of 5 X lO'^My, from Daddi et al. (2007a) 
and Elbaz et al. (2007). The dashed lines show the measured l-cr range (about 0.2 dex at z = 0.1 and about 0.3 dex at z = 1 and 1.9 and 0.36 dex at z = 3.8). We 
recall that stellar masses and SFRs in this figure are shown for the Chabrier (2003) IMF. 



photometry is ^ 2.3 x 10 M©, expressed for a Chabrier 
(2003) IMF (see, e.g., Maraston et al. 2006 for more details 
on the range of models and assumptions implied). The com- 
bined molecular gas and stellar masses, only roughly half of 
which presumably fall within r^, add up to roughly 70-80% of 
the estimated dynamical mass, close enough within the large 
uncertainties of such estimates. This implies that the stellar 
mass is not highly overestimated for the choice of a Chabrier 
(2003) IMF and the case of a top heavy IMF is disfavored by 
the data. 

For GN20.2a, we infer with a similar method a stellar mass 
of -0.5 X 10^ 'Mq. 

4.4. On the possibility oflensing 

The intrinsic SFRs and luminosities of these sources might 
be lower if there is amplification by gravitational lensing (see 
Paciga & Scott 2008). Lensing by foreground structure could 
be an alternative way to interpret the fact that colder SEDs, 
typical of lower luminosity galaxies in the local Universe, 
seem to reproduce better the colors of the GN20 and GN20.2 
galaxies. In that case, also the luminosity would match well to 
those of colder galaxies, without the need to advocate evolu- 
tion in the temperature-luminosity relation of starburst galax- 
ies. Some evidence against lensing is however provided by 
the FIR to CO properties of GN20 and GN20.2 (Fig. [IB. If 
lensed by a factor of 10 or more, the de-lensed properties 
would place them in a region of the diagram corresponding 
to a very high star formation efficiency, not commonly ob- 
served for LiR 1O'^L0 sources at high and low redshifts. 
Following the analysis discussed in Tacconi et al. (2008), the 
agreement between the dynamical mass estimate and its stel- 
lar and gas mass estimates also disfavors the possibility that 



GN20 is lensed by a very large factor, given that emission 
line FWHMs should be independent on magnification. Qual- 
itatively, the same picture also holds for GN20.2a, for which 
the data suggest a CO[4-3] hne FWHM- 1000 km s"' (albeit 
with a large error), not uncommon amongst SMGs (Greve et 
al. 2005). 

4.5. The relation between stellar mass and SFR at z ^ 4 

A tight correlation exists between the SFR and stellar mass 
of star forming galaxies from z = to z — 2, after exclud- 
ing the locus of quiescent/passive galaxies (Elbaz et al. 2007; 
Daddi et al. 2007a; Noeske et al. 2007). The implication of 
this correlation remains subject to debate but suggests that, on 
average, the processes regulating the star formation activity 
of a galaxy are the same over a large range of stellar masses. 
Outliers do exist at both extremes in this correlation. On the 
low SFR side, the cloud of red-dead galaxies have SFRs lower 
than expected by several orders of magnitudes. On the high 
SFR side, ULIRGs and z — 2-3 SMGs exhibit approximately 
one order of magnitude larger SFRs for their stellar mass than 
the average galaxy of equivalent stellar mass (Daddi et al. 
2007a; Elbaz et al. 2007). Here, following Daddi et al 2007a, 
we use the Tacconi et al. 2008 estimate of SSFR for z — 2-3 
SMGs that are based on dynamical masses rather than stellar 
masses which should be a more robust derivation of the total 
galaxy masses, and we converted the average Lir ~ IO^^ 'L© 
of the sample into SFR - 1200 M© yr"'. Takagi et al. (2008) 
confirm higher SSFRs in SMGs, respect to normal galaxies 
of the same masses, using stellar mass estimates, although 
finding SSFRs a factor of 3 higher on average and a few pos- 
sibly overlapping galaxies (see also Dannerbauer et al. 2006 
for overlapping objects between the BzK and SMG samples). 
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Their SFRs are based on the radio and assume a radio-IR cor- 
relation with a different normalization respect to the local one 
by a factor of about 2. 

On the low side, the red-dead galaxies appear to fall short of 
SFR because of a lack of molecular gas fuel, possibly due to 
negative feedback, the excess SFR of z ^ 2-3 SMGs remains 
to be understood. A possible explanation is external trigger- 
ing, e.g. by major mergers (Tacconi et al. 2006; 2008) as it is 
seen for local ULlRGs. 

When compared to their closer z ^ 2-3 siblings, z ^ 4 
SMGs appear to be forming stars with an equivalently high 
SSFR. However, the stellar mass-SFR relation shows a con- 
tinuous increase of SSFR with increasing redshift, which is 
expected due to the higher gas mass fractions at higher red- 
shifts. The question therefore arises as to whether an SMG 
could represent the typical star forming galaxy at z ^ 4, con- 
trary to their z ^ 2 siblings which are atypical compared to 
the average galaxy at their cosmic epoch, having much larger 
SSFR, see Fig. [H] 

In order to address this question, one needs to identify the 
locus of average star forming galaxies at z ^ 4 and then com- 
pare it to the locus of z 4 SMGs. We consider the sub- 
sample of B-band dropout Lyman break galaxies in GOODS- 
N with IRAC magnitudes at 5.8 and 8.0/im brighter than 24. 
The requirement of an IRAC detection is necessary to reli- 
ably estimate the stellar masses. The use of B-band dropouts 
ensures that we are concentrating on star forming galaxies, 
with a negligible contribution from passive galaxies. Within 
this limit, about 10% of the Z?-band dropouts to zab < 27 are 
selected, or 77 sources. We fitted the ACS to IRAC SEDs 
of these B-band dropouts using Maraston (2005) models with 
constant star formation, allowing the redshift to vary freely 
within 3.1 < z < 4.4 (the 95% confidence region expected 
for the redshift distribution of B-band dropouts) and allow- 
ing also for dust reddening following a Calzetti (2000) law. 
In this way, we estimated stellar masses and (basically UV- 
driven) SFRs for each of the B-band dropout Lyman break 
galaxy. Typical errors from the fit amount to factors of about 
2 for both quantities. From this we derive an average stellar 
mass ofl.5xlO'°Meand an average SFR of 4OM0 yr We 
performed stacking at 850/im using the SCUBA super-map of 
Pope et al. (2005) and a technique identical to what was done 
in Daddi et al. (2005b; 2007a). We do not detect the galaxies, 
deriving a 3(7 flux density limit of 0.54mJy, corresponding to 
57 Mq yr^^ From radio stacking, again performed with a 
similar technique to Daddi et al. (2005b; 2007a), we get a 
tentative detection with a peak flux of 1.7 ±0.6/iJy. The 3(7 
radio upper limit corresponds to a SFR of 54M0 yr^', consis- 
tent with the limit from submm stacking. Overall, the stack- 
ing limits appear to be consistent with the SFRs derived from 
SED fitting in these galaxies. This is similar to what found at 
z 2 where analogous estimates based on the UV luminos- 
ity of galaxies appear to agree well on average with mid- and 
far-lR estimates (Daddi et al. 2007a). 

The locus of typical, massive z ^ 4 B-band dropouts does 
not support a continuously increasing SSFR with redshift, 
suggesting instead a plateau of the SSFR for z ;^ 2. The aver- 
age SFR of 40 Mft yr^^ places these 'typical' z '--^ 4 galaxies 
fairly close to the z ~ 2 correlation. 

The B-band dropout Lyman break galaxies exhibit an aver- 
age stellar mass lower than the z ~ 4 SMGs. This implies that 
not only do z ~ 4 SMGs present similar SSFRs to their z 2 
siblings, but they are also found to be outliers with respect 



to the average star-forming galaxy at their epoch. Note that 
the similarity of z = 2-3 and z = 4 SMGs is also supported 
by their equivalently large FIR to CO ratios, corresponding 
to high star formation efficiencies and rapid gas consumption 
timescales. 

We hypothesize that this favors a triggering mechanism for 
the activity of the z ~ 4 SMGs, such as major mergers. In- 
deed, GN20.2a has a very nearby (0.5" or 3.5 proper kpc at 
z — 4.05) B-band dropout companion (Fig. |5]l. GN20 has a 
clumpy morphology that could be reminiscent of a merger (al- 
though, see Bournaud et al. 2008). The 0.5" offset between 
the optical light peak and submm/radio and CO might also 
point to a complex geometrical situation often found among 
mergers (e.g., the Antennae; Wang et al. 2004). The high gas 
densities of a merging event could also explain why, contrary 
to what is found for typical massive galaxies at 1 .4 < z < 2.5 
(Daddi et al. 2007a; Dannerbauer et al. 2006), the emitting 
region is completely opaque in the UV, similar to z 2-3 
SMGs. Indeed, we find that the radio-lR based SFR estimates 
of the z 4 SMGs are a factor of > 10 times larger than the 
UV-based ones. 

5. A z = 4.05 PROTO CLUSTER OF GALAXIES IN GOODS-N 
CENTERED AROUND GN20 AND GN20.2 

Down to the limits of the GOODS-N observations of P06, 
the surface density of SMGs is w 1000 deg^^, and thus the 
fact that GN20 and GN20.2 are only 24" from each other, 
with a projected spatial distance of 169 physical kpc, and both 
are at z ~ 4.05, implies that these two SMGS are spatially 
clustered (see also Blain et al. 2004). This suggests that these 
SMGs might also be part of some type of larger, high redshift 
structure in GOODS-N. 

To investigate further this possibility, we have analyzed 
the distribution of spectroscopically confirmed Lyman-break 
galaxies at z 4 (B-band dropouts) in GOODS-N (Fig. [131 
Stern et al., in preparation). We compare the observed redshift 
distribution to the expected distribution function derived from 
simulations (Giavalisco et al. 2004b; Lee et al. 2006). There 
are 12 redshifts in the Az = 0.04 bin that includes GN20 
(z = 4.055); these galaxies define a clear spike in the red- 
shift distribution. There currently are 58 high quality spectro- 
scopic redshifts of B-band dropouts in the GOODS-N sample 
distributed over the range 3.5 < z < 4.5, 12 of which belong 
to the spike. This is a very strong concentration of galaxies, 
comparable, for example, to that discovered at z '--^ 3.09 by 
Steidel et al. (1998), and it seems unlikely that it can be the 
result of random (Poisson) sampling. We have generated a 
large (10^) ensemble of random realizations of GOODS B- 
band dropouts of the same size as the observed one using the 
expected redshift distribution function, and have counted how 
often a spike as large as the observed one is realized. In prin- 
ciple, one should also take into account the spatial clustering 
of B-band dropouts in this simulation; in practice, however, 
the spatial correlation length, ^ A comoving Mpc, of these 
sources (Lee et al. 2006) is 7 x smaller than the radial distance 
covered by the redshift bin that we have adopted, Ar — 29 co- 
moving Mpc, implying that over these scales the spatial dis- 
tribution of the galaxies can be approximated as Poisson. 

We found no realization in which the number of galaxies 
in the redshift bin of GN20 equals or exceeds that of the 
observed spike, meaning that the probability to get a simi- 
lar overdensity by chance is less than 10^^. The simulations 
show that the average number of galaxies in the 0.04 redshift 
bin is 1.31 ± 1.18, implying that the observed spike corre- 
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Fig. 13. — Left panel: spectroscopic redshift distribution of GOODS-N B-band dropout Lyman break galaxies galaxies, in bins of 0.04 (solid histogram). A 
strong spike is seen, coincident with the redshift of GN20 and GN20.2a (vertical dashed line). The dotted histogram shows the expected redshift distribution 
of B-band dropouts, based on the simulations of Giavalisco et al. (2004b) and Lee et al. (2006). Right panel: spatial distribution in the northern comer of the 
GOODS-N field (about 3.5' x 3.5') of zab < 27 galaxies (small points), and B-band dropouts (crosses). The filled symbols are GN20, GN20.2a and GN20.2b. 
The large circle shows a 25" (180 kpc) radius area centered on GN20. The blue filled symbols are two spectroscopically confirmed galaxies at z = 4.058 (closer; 
this is BD29079, also shown in Fig.|5) and z = 4.047 (further). 



sponds to an overdensity of a factor of 8.4 and significant at 
the 8.2(7 level. These redshift concentrations are frequently 
found among Lyman-break galaxies at z ^ 3. For example, 
similar concentrations are found in every survey field for U - 
band dropouts by Steidel et al. (2003). Vanzella et al. (2007) 
also report evidence of concentrations in the redshift distribu- 
tion of B-band dropouts in the GOODS-South field. 

To gain additional information about the nature of the GN20 
structure, especially its transverse size, we have studied the 
distribution of the angular positions of B-band dropouts in 
the field and around the GN20 complex. Of the 765 Z?-band 
dropouts to zab < 27 in GOODS-N, counting all sources re- 
gardless of a spectroscopic identification, 14 are within 25" 
from GN20 (2.5 Mpc^ comoving; see Fig. fT4l l. This is an 
overdensity of a factor of 5.4, the strongest concentration 
of B-band dropout Lyman break galaxies being present in 
GOODS-N. Given that at these angular scales the angular cor- 
relation function of B-band dropouts is non negligible (Lee et 
al. 2006), we account for the variance due to clustering to de- 
rive an expected rm.s fluctuation in the counts in 25" radius 
apertures of about 1.8 galaxies. Therefore, the excess counts 
that we see in the 25" radius region surrounding GN20 corre- 
spond to a 5.8(7 fluctuation. Taken together with the analysis 
of the redshift spike, this leads us to conclude that the GN20- 
GN20.2 SMG pair is part of a very strong overdensity with a 
size consistent with a cluster of galaxies, most likely a proto- 
cluster given the redshift. 

Additional evidence for the presence of a strong overden- 
sity in this structure comes from Mancini et al. (2008), 
who use photometric redshifts to identify galaxy candidates at 
z > 3.5 from an IRAC sample limited to m4 < 23. They 
find a total of 54 candidates in their full GOODS-N IRAC 
catalogs, which covers 180 arcmin^. Of these massive galax- 
ies, four are within 25" from GN20 (GN20 itself, GN20.2a, 
GN20.2b and BD29079; all are shown in Fig. H). This cor- 
responds to an overdensity of a factor of 18 for this IRAC- 



selected sample, with a chance probability of ^ 10^^ assum- 
ing no spatial clustering (Poisson statistics; 0.16 galaxies ex- 
pected on average and three observed in addition to GN20). 

The possibility that we have uncovered a proto-cluster is re- 
inforced by the fact that a number of massive galaxies appear 
to be part of this structure. With an estimated total baryonic 
mass of ~ 2 - 3 X lO'^M©, GN20 is likely one of the most 
massive galaxies of this structure. The presence of over a 
dozen additional B-band dropouts in this 2 Mpc^ spatial con- 
centration implies a total baryonic mass reaching ^ lO'^M© 
for the overdensity. If we assume the typical mass-to-light ra- 
tio of local massive clusters of ^ 50 (Lin, Mohr & Stanford 
2003), this suggests a lower limit to the total (dark) mass of 
the proto-cluster of about 5 x IO^^Mq. 

These lines of evidence all support the identification of a 
large structure at redshift z ^ 4. While its transverse size 
appears to be around 2 Mpc (comoving) in diameter, its ra- 
dial size remains basically unconstrained given that the er- 
rors on the redshift measurements are already comparable to 
the transverse size, not even considering peculiar motions 
along the radial direction. In addition, a significant over- 
density at z = 4.05 seems to be present extending over the 
whole GOODS-N field. This suggests that we are witnessing 
a proto-cluster of galaxies which, considering the incomplete- 
ness in all the samples we have been discussing, encompasses 
a total mass of « 

It is interesting to compare the GN20 environment with 
other structures previously discovered at high redshift around 
powerful steep-spectrum radio galaxies (Kurk et al. 2000; 
Venemans et al. 2002; 2007; Miley & de Breuck 2008 and 
references therein). The most direct parallel is with the proto- 
cluster around TN J1338-1942 (De Breuck et al. 1999; Ven- 
emans et al. 2002) at z 4.1. A study of the overdensity 
of B-band dropout galaxies around this source has been per- 
formed by Miley et al. (2004), reporting factors of 2.5 — 5 ex- 
cess of Lyman break galaxies (significant at the 3 — 5(7) level. 



Two bright submillimeter galaxies in a z = 4.05 proto-cluster in GOODS -North 



17 



N 



BzK-21000 



GN20' 



BD29079 



.4 : 



GN20.2a 



J 



■ / 

GN20.2b 



/ 



Fig. lA.— HST+ACS color image of the field around GN20, derived from F435W (blue), F775W (green) and F850LP (red). The image size is 50", or 1.77 
comoving Mpc at z = 4.05. Close-ups of GN20, GN20.2a, GN20.2b and the luminous "B-drop" are shown in Fig. [3] The BzK-21000 galaxy at the top of the 
image was the original target of the PdBI observations (Daddi et al. 2008). Additional dropouts in the field can be easily recognized from their yellow colors 
(e.g., see Fig. 12). 



depending on the spatial scale. This is fairly similar to what is 
found here for the overdensity around GN20. The field around 
TN J1338-1942 was also observed at 1.2 mm by de Breuck et 
al. (2004), finding an overdensity at mm wavelengths as well, 
although without spectroscopic confirmations so far 

Also, it is interesting to note that, contrary to the structure 
described by Chapman et al. (2008a) at z = 1.99 where 5 
SMGs are found but without a strong overdensity of UV se- 
lected galaxies, in this case the two SMGs that we have found 
at z = 4.05 (and more SMGs might well be at z = 4.05 in the 
GOODS-N field, see next Sect.|6]l do appear to coincide with 
a strong overdensity of Lyman break galaxies as well. 

6. RADIO-IR PHOTOMETRIC REDSHIFTS FOR SMGS 

An important result of this paper is the success in identi- 
fying the z 4 redshifts for the GN20 and GN20.2 galaxies 
based only on the flux measurements at 24/im, 850/im, and 
20cm. The ratio of 20cm to 850//m flux density is, alone, a 
powerful redshift indicator for dusty submm galaxies (Carilli 



& Yun 1999; 2000)'"*, although Hmited in its application by 
its strong dependence on SEE) temperature. Clearly, colder 
galaxies (P06) are intrinsically brighter at 850/im at fixed ra- 
dio flux density and redshift (e.g., P06). In our approach, 
however, we have exploited the 24/im measurement that lies 
on the opposite side of the SED peak. It turns out that for the 
variations in the shape of local SEDs (as codified in the CEOl 
library) the 850/im to 24/im flux density ratio is fairly con- 
stant at fixed redshift, compared to the redshift excursion of 
its value (Fig.fTSll — e.g., this ratio changes by factors of ^ 2 
at fixed redshift, while changing the bolometric luminosity 
LjR by a factor of 10. However, at fixed Ljr this ratio changes 
by almost four orders of magnitude when going from z — 0.5 
to z = 5 (Fig. [TSl l. Note that we have not re-normaUzed the 

We applied this redshift indicator to the counterparts of the two SMGs 
presented here. For GN20, we obtain z = 3.74. For GN 20.2a, we obtain 
z= 1.61 and for GN 20.2b we obtain z = 4.04. The low photometric redshift 
obtained for GN 20.2a indicates that the radio emission must be contaminated 
by an AGN, which substantially affects the Carilli & Yun indicator. 
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Redshift 

Fig. 15. — The redshift dependences of several flux density ratios and for 
the 1.4 GHz radio flux density expected for CEOl models with Lir = 10"L,.j 
(dashed line), Lir = lO'-L... (solid line) and Ljr = lO'^^Lv, (dotted line). The 
filled symbols show the values for GN20, which would be detectable up to 
z ~ 8 to the 5-cr limits of 20/^Jy in our 1 .4 GHz map in GOODS-N (Morrison 
et al., in preparation). 



CEOl template models when fitting observed data from the 
GN20 and GN20.2 galaxies (Fig.O bottom panel). This in- 
formation is a powerful complement to the flux density ratios, 
as it helps pinpoint the correct intrinsic source luminosities, 
thus helping to solve dust temperature degeneracies, assuming 
that the local correlations between temperature and luminos- 
ity measured by CEOl and the radio-FlR correlation are not 
severely altered in the distant universe. As a result, the photo- 
metric redshift estimates based on IR and radio appear to be 
robust against SED temperature/shape variations, as shown by 
the case of GN20 and its possibly colder SED (Fig.|8]and [TOl i. 

Encouraged by the apparent success of this technique, we 
explored in more detail the possibility of deriving accurate 
radio-lR photometric redshifts using mid-lR, submm, and ra- 
dio observations. We applied this technique to the sample 
of well-studied SMGs in GOODS-N (P06) and the Lockman 
Hole (Scott et al. 2002; Ivison et al. 2002; Egami et al. 
2004), considering only the most likely counterparts. Our 
radio-lR photometric redshifts are determined as described in 
Sect. 13.1.21 for GN20. In particular, in order to effectively 
use the 850/xm information in the fit, we limit the S/N of the 
24/xm and 20cm flux density to a maximum of 10, artificially 
increasing the relative error bars if the S/N is higher Table |4] 
summarizes the radio-lR photometric redshift constraints. 

Results of the comparison are in Fig. [16] In the left panel, 
we show SMGs with counterparts having a measured spectro- 
scopic redshifts. The agreement between radio-lR photomet- 
ric redshifts and spectroscopic redshifts are in general fairly 
good, with a semi interquartile range of Az/(1 +z) — 0.12. 
This figure includes also the Capak et al. (2008) source (as- 
suming a ratio of 2.8 between the 850/im and 1.2mm flux 
densities), but the result does not change substantially if this 
one is excluded. 
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Table 4. Radio-IR photometric redshifts of SMGs 



ID 


ZORIG 


zType 


radio— IR 
^best 


radio-IR radio-IR 2 
^min MTiax X 


GNOl 


2.415 


2 
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19^ 
i .Z J 


Z. Jj 


1 70 

i. ly 


GN()6 


2.00 


2 


2.30 


z. iU 


T /in 
Z.4U 


n /I Q 


GN()7a 




_! 


3.65 


J. Jj 


A 1 n 


nz 


GN()7b 


1.998 


2 


1.95 




T on 
Z.ZU 


ifS. Ij 


GN()9 


2.9 





13.42 


z.yu 


1 ^ ^n 
1 J. jU 


1 Q 


GNIO 


2.2 





5 43 


3.75 


1 1.90 


3.02 


GNl 1 


2.3 





2 45 


2.20 


2.60 


17.89 


GN12 


3.10 


I 


3 20 


2.90 


3.50 


n AO 


GN13 


475 


2 


1 30 


1.05 




A 1 c 

4. iJ 


GN15 


2 74 


I 


2 15 


1 .90 


2.35 


23.27 


GN16 


1.68 


I 


2 50 


2.40 


2.60 


12.74 


GN17 


1 73 


2 


1 35 


1.15 


l.uD 


5.98 


GNl 8 


2 8 


Q 


2 50 


2.25 


2.70 


1.4 J 


GNl 9a 


2 484 


3 


3 15 


2.70 


6.00 


n zc\ 


GNl 9b 


2 484 


3 


1 35 


0.05 


2. 10 


5.23 


GN20 


4 055 


3 


4 20 


3.45 


4.80 


17.59 




4 051 


3 


3 70 


i.jj 


A 

4.Uj 


Zl .Z / 


GN21 


2 8 


Q 




2.60 


3.60 


n oc 
O.ZD 


GN22 


2 509 


2 


2 95 


2.60 


3.50 


1.21 


GN23 


2.6 





2.95 


Z.J J 


1 1^ 

D.JJ 


1/10 

i.4o 


GN24 


2.91 


1 


2.15 


U. iU 


z. /u 


J. io 


GN25 


1 013 


2 


1 15 


(\ 1^ 

U.JJ 


1 /in 


1 /I Q 


GN26 


1.223 


3 


1.82 


1 

1 . JU 


i.Vj 


1 1 Q/1 


GN28 


1.020 


2 


2.95 


1 

Z.J J 


i4.oU 


1 n 1Q 
iU.Jo 


GN30 


1.355 


2 


2.15 


1.45 


2.45 


0.37 


GN31 


0.935 


2 


0.90 


0.05 


1.30 


8.86 


GN32 


1.9 





1 75 


0.05 


2.35 


4.84 


GN37 


3.190 


2 


2.90 


0.10 


15.50 


1.89 


LE850.1 


2.60 


1 


2.30 


2.10 


2.80 


10.65 


LE850.4 


2.60 


1 


6.18 


3.00 


15.50 


0.02 


LE850.7 


1.80 


1 


2.20 


1.95 


2.45 


1.18 


LE850.8b 


3.00 


1 


2.10 


1.35 


2.55 


2.96 


LE850.8a 


0.974 


2 


0.93 


0.05 


6.20 


14.83 


LE850.14b 


2.50 


1 


2.35 


2.05 


3.05 


0.56 


LE850.14a 


2.380 


2 


2.55 


1.45 


4.35 


0.01 


LE850.18 


2.690 


2 


2.45 


1.95 


3.15 


0.28 


LE850.35 


3.00 


1 


2.40 


1.85 


2.80 


1.11 


Capak 


4.547 


3 


3.97 


3.00 


5.20 


0.22 


NOTE. — GN 


galaxies 


are from P()6. LE850 gal 


dxies are 


from Egami 


et al. (2004) 



The last object is from Capak et al. (2008). 'zType' indicates the origin of stellar pho- 
tometric redshifts (cqrig): for 5/? (7ze/-l-IR AC/MIPS photometric redshifts (P06), 1 for 
ordinary photometric redshifts from stellar emission, 2 for spectroscopic redshifts (usu- 
ally optical/UV redshifts or Spitzer IRS redshifts derived mainly from detection of PAH 
features), 3 for CO redshifts. The CO redshifts for GN19a and GN19b are from Tac- 
coni et al. (2008), the one for GN26 is from Prayer et al. (2008), those for GN20 and 
GN20.2a are from this paper. Some of the P06 photometric redshifts listed in P06 were 
slightly off; here we list the correct values. New spectroscopic redshifts for three P06 
galaxies were measured with IRS in Pope et al. (2008). Ranges in r-raiiio-iR gj-g given 
at the 99% confidence level (Ax^ — 6.63) following Avni et al. (1976) for a single 
interesting parameter in the fit (e.g., redshift). For the cases of double components we 
assigned the total submm flux to each component while computing the photometric red- 
shift. The photometric redshifts listed here for the GN20 and GN20.2a galaxies did not 
use the 3.3mm continuum measurements, for consistency with the rest of the galaxies. 



We identify three galaxies as being suspect counterparts to 
the submm emission based on their 99% confidence error bars 
falHng outside the median photometric redshift relationship 
(Fig. 15). These three galaxies are from the P06 sample, and 
they are all characterized by being at 6'' separation or more 
from the SCUBA position. Given that P06 finds a median 
separation of ^ ?>" for secure counterparts, with 80% being 
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Fig. 16. — Left panel: radio-IR photometric redshifts for SCUBA counterparts with known spectroscopic redshifts from P06 for GOODS-N and Lehmann et 
al. (2001), Ivison et al. (2005), Chapman et al.(2005) for Lockman Hole. Empty triangles with error bars (99% levels) are GN13, GN04.2 and GN28 from P06. 
We suggest that these sources are not the correct counterparts to the submm emission based on the discrepant radio-IR photometric and spectroscopic redshifts. 
Right panel: the same comparison, but with counterparts having only photometric redshifts determined from the properties of the stellar SED, from P06 and 
Egami et al. (2004). Tiiangles with en'or bars are sources which are consistent with a z > 4 solution of the radio-IR photometric redshifts. 



at < 4.5" , the large offsets alone suggest that the counter- 
part identification might have actually failed for these objects. 
Only three out of 43 targets are identified as such, and two of 
these were already listed as tentative counterparts, confirming 
the high quality of counterpart identification in the original 
work. We discuss these three objects here in some detail: 

• GN13: the counterpart in P06 is a galaxy at Zspec — 

0.475 

lying at 7" from the submm position. The identification was 
judged to be secure, on the basis of the low probability of a 
chance association between a radio and mid-lR bright source 
relative to the submm position. Its 850/im flux density of 
1.9 mJy, despite being quite faint intrinsically, is about six 
times brighter than expected on the basis of both the 24/im 
and 20cm radio flux densities for its redshift, assuming the 
CEOl templates. We inspected the other galaxies present in 
the submm beam, and found no other convincing counterpart. 
This SMG could be at quite high redshift. 

• GN04.2: the counterpart listed by P06 is at Zspec = 0.851 and 
5.7" from the submm position. The 24/im and radio flux den- 
sities are both too faint to be consistent with the 850/im flux 
density of 2.7 mJy at this redshift, assuming the CEOl tem- 
plates. Furthermore, this source is no longer detected in the 
latest and deepest 20cm reduction (Morrison et al., in prepa- 
ration); the P06 results were based on a 3(t detection at 20cm. 
No other 24/im detected galaxy is present in the beam, which 
may imply that the real counterpart is at fairly high redshift. 

• GN28: the counterpart listed by P06 is at Zspec = 1.020 and 
at 5.8" from the submm position. There is a strong incon- 
sistency between the faint 24/im flux density of this galaxy 
(20.5/iJy) and its 850/im emission (1.7 mJy). The 850/im to 
24/im ratio is two orders of magnitude larger than expected 
for starforming galaxies at z = 1 (Fig.fTSll: and would be even 
more discrepant for an AGN SED. The 20cm radio position 

" Dannerbauer et al. (2004) also finds a median separation of 3" for the 
radio counterparts of MAMBO sources, with a similar percentage of radio 
counterparts being separated by less than 4.5". 



is close to the radio lobes of a bright radio galaxy at a simi- 
lar redshift, and thus may be severely contaminated. No other 
24/im detected galaxy is present in the beam, which may also 
imply that the real counterpart is at fairly high redshift. 

When excluding these three likely incorrect counterparts, 
we find (t(Az/(1 +z)) = 0.15 and a semi-interquartile range 
of Az/(l+z)=0.08. 

The right panel of Fig. [16] shows the comparison of the 
radio-IR photometric redshifts with those counterparts of 
SMGs from the same works that have only photometric red- 
shift estimates based largely on the stellar emission (plus the 
IRAC to MIPS colors in a few cases from P06). Again, the 
agreement is quite good, with a similar cr(Az/(l +z)) = 0.14. 
We find only two cases of widely discrepant estimates, dis- 
cussed here in some detail: 

• GNIO: P06 assign a z = 2.2 photometric redshift, while 
we find z > 3.7 at the 99% confidence level, in good agree- 
ment with Wang et al. (2007; 2008) and Dannerbauer et al. 
(2008). The discrepancy is likely due to the photometric red- 
shift of P06 being too low, while the counterpart identification 
appears correct and agrees with Wang et al. (2007; 2008) and 
Dannerbauer et al. (2008). 

• GN16: The P06 identification appears correct, and the dis- 
crepancy is due to the low photometric redshift (Zphot = 1 -68) 
estimate of P06. 

We have also applied the radio-IR photometric redshift 
technique to the CO-detected SMG galaxy reported by Capak 
et al. (2008) and Schinnerer et al. (2008) to be at z = 4.547. 
The resulting photometric redshift is z = 4.0 and the 99% con- 
fidence range is 3.0 < z < 5.2, supporting the z = 4.5 red- 
shift identification of this SMG despite the large offset ( ^ 1"; 
}^ 10 kpc proper) between its Lya emitting optical counter- 
part and the radio emission. 

In conclusion, we find that radio-IR photometric redshifts 
can be used as a powerful and independent complement 
to photometric redshifts derived from stellar emission, and 
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Table 5. Properties of the z > 3.5 starburst candidates selected using 
their high radio 1 .4 ghz to flux density ratios 



ID RA(J2000) DEC(J2000) Si.4GHz 524^1, Ss.OMm 5'4.5Mn 

(/iJy) (p}y) (^l3y) (^l3y) 



VLA512 12:36:58.49 62:09:31.8 47.3 ± 11.5 36.8 ± 5.0 9.7±1.0 4.0±0.4 
VLA812 12:35:53.24 62:13:37.5 45.6 ±5.1 33.0 ± 6.2 23.1 ±2.3 14.6±1.5 
VLA966 12:37:01.52 62:20:24.6 64.2 ± 5.3 49.7 ± 3.6 22.9±2.3 14.2±1.4 



could greatly improve the SMG counterpart identification pro- 
cess. Some differences and evolutionary trends are seen when 
studying the SEDs of SMGs in detail (e.g.. Pope et al. 2006; 
2008). However, the fact that direct fitting of locally-derived 
CEOl templates yields reasonable photometric redshifts with- 
out any significant bias with respect to redshift suggests that 
ULIRGs SEDs have not dramatically evolved from z = to 

7. THE SELECTION AND SPACE DENSITY OF THE HIGHEST 
REDSHIFT STARBURST GALAXIES 

We emphasize that not all radio faint (e.g., 5i.4ghz < 
30/iJy) counterparts to SMGs are likely to be at z > 4. For 
example, four SMGs (GN03, GNU, GNU (HDF 850-1), 
GN18) have even fainter radio flux densities than GNIO and 
have similar properties to GNIO at other wavelengths, e.g. 
they are undetected in ACS imaging and have red IRAC 
SEDs. Dannerbauer et al. (2008) suggest they could lie at 
z ^ 4. In contrast, the radio-IR photometric redshifts suggest 
that GN03, GNU and GN18 are at 2 < z < 3, which is still 
consistent with their IRAC SEDs. This is because these galax- 
ies have quite bright 24/^m flux densities, much larger than 
the radio ones. Besides GN20, GN20.2a and GN20.2b, only 
four additional objects from P06 are identified by the photo- 
metric redshift technique to possibly be at z ^ 4. With vari- 
ous degrees of confidence, they are GNIO, a very solid z ^ 4 
source also identified by Wang et al. (2007; 2008) and Dan- 
nerbauer et al. (2008), and GN13, GN28 and GN04.2 (where 
reliable counterparts are, however lacking, thus making these 
less convincing candidates). Another well known z 4 candi- 
date from GOODS-N is the galaxy HDF850.1 (Dunlop et al. 
2004), listed as GN14 in P06 for which we have not obtained 
a radio-IR photometric redshift due to the lack of reliable ra- 
dio and 24/im data. This would make a total of eight plau- 
sible z > 3.5 SMGs in GOODS-N in a total ai-ea of roughly 
50—100 arcmin^ covered with SCUBA (with a highly inho- 
mogeneous coverage; see, e.g.. Pope et al. 2005), or roughly 
200 — 300 galaxies deg^^. Using the volume in 3.5 < z < 6, 
this corresponds to a space density of « 10^^ Mpc^^, which 
compares fairly well with the expectations for the very high 
redshift formation of the oldest early type galaxies seen at 
z ^ 1.5 — 2.5 (see Sect. 1). With typical SFRs in excess of 
thousand Mq yr^^, such a population contributes an SFRD of 
order of 2 x 10^^ Mq yr^' Mpc^^, already comparable to the 
global contribution of Lyman break galaxies at these redshifts 
(Giavalisco et al. 2004b). 

7.1. A radio/24 fim selection criterion for high-z starburst 

galaxies 

For the most reliable z^4 starbursts (GNIO, GN20, 
GN20.2a) we notice that all of these have very high 20cm 
to 24/im flux density ratios, greater than about 1 (Fig. [TsT l. 



GN20.2b itself satisfies such a criterion and might also be 
at z = 4.05. This suggests that a Sia ghz ^ S24fj,m criterion 
could be useful at singling out the very high redshift tail of 
SMGs. Applying this criterion to the P06 sample we re- 
cover the SMGs Hsted above, plus GN12, GN16 and GN19a, 
all sources at 2.5 < z 3.2, a distribution clearly skewed to 
higher redshifts than typical SMGs. The z = 4.547 galaxy 
from Capak et al. (2008) also satisfies this criterion. 

In principle, this simple criterion could be applied regard- 
less of a submm detection, and could avoid the inherent bias 
toward intrinsically cold sources. Of course, care should be 
taken and a careful analysis of all multiwavelength informa- 
tion should be performed on targets selected in this way, as 
radio emission from radio-loud AGN would bias the radio 
to 24/im flux density ratio. Mid-IR as well can be boosted 
by AGNs, which could make genuine z > 4 objects fail this 
criterion. More importantly, given the relatively small dy- 
namic range expected over redshift for the 20cm to 24/im 
flux density ratios in galaxies (Fig. [TSl l. scatter due to SED 
variations can imply a substantial contamination from lower 
redshift sources. 

We have performed a preliminary test of this criterion with 
a sample of radio selected galaxies with 5'i.4 ghz > 30/iJy in 
GOODS-N (Morrison et al., in preparation), applying a sim- 
ple but very conservative approach. By requiring Si a ghz > 
0.8 X S24fj.m, 63 radio sources are retained. Only half of these 
(34) are undetected in the ACS F435W band (S/N < 2), a 
necessary condition for being at z^4-. Of these, 16 are also 
consistent with a red IRAC SED, with a peak beyond Ifim, 
as expected for z^4'^. We further require relatively blue 
5p/fze r+MIPS to IRAC flux density ratios: 524^m/'S4.5^m < 13 
and 5'24^m/ •S'g o/jm < 5, as observed in the bona-fide z ^ 4 sam- 
ple (GN20, GN20.2a, GN20.2b, GNIO), leaving us with 11 
galaxies. This excludes from the sample the z — 4.424 radio 
source from Waddington et al. (1999), where the near-IR to 
mid-IR spectral range is likely dominated by a powerful AGN. 
Finally, we exclude all sources detected in the X-ray catalog 
of Alexander et al. (2003) which, if at z ^ 4, must be AGN 
dominated. This leaves us with seven sources, including the 
three reliable z > 4 galaxies already known (GN20, GN20.2a, 
GNIO) plus GN20.2b. We consider the properties of the three 
new candidates in detail: 

• VLA 512: this object is undetected in a relatively shallow 
region of the SCUBA supermap. However, some signal is 
found at its position with S'l.imm = 3.2mJy (S/N ^ 3. 4-) in an 
AzTEC map of GOODS-N (G. Wilson, private communica- 
tion). This source has no counterpart at any ACS band, and is 

The Capak et al. (2008) galaxy does not satisfy this requirement. The 
authors suggest that line emission contaminate the IRAC photometry for this 
source. A possible alternative is that the IRAC flux detected is actually from 
a foreground, lower redshift galaxy. This would imply that the SMG is very 
faint even at the IRAC bands, perhaps similar to GNIO. 
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a good high redshift starburst candidate. 

• VLA 812: this galaxy is very close to the edge of the 
SCUBA supermap, where the noise is fairly high. However, 
this source is listed in Chapman et al. (2005) as a SCUBA 
detection in photometry mode (^gso^m = 8.8 ± 2.1mJy), with 
a redshift of z = 2.098 (the source is also detected by Per- 
era et al. 2008 and Greve et al. 2008 at 1.1mm and 1.2mm 
with MAMBO and AzTEC, respectively). We find that a rel- 
atively bright z — 2.098 galaxy is present in the area, but 2" 
away from the IRAC and radio position, which corresponds 
instead to a faint galaxy with ACS photometry of /ab = 27.6 
and ZAB = 27.3, undetected at B- and v-band. We suggest 
that the Chapman et al. (2005) redshift identification does not 
correspond to the actual counterpart of the radio and submm 
emission. This radio source is a reliable candidate for a very 
high redshift starburst (z ^ 4). 

• VLA 966: we place a 3cr upper limit of S%5Q^rn < 8.1mJy 
from the SCUBA supermap at the position of this source. Its 
optical counterpart is a B-band dropout with z = 26.6 and red 
IRAC SED, confirming that this is a reliable z ~ 4 candidate. 
The VLA radio flux density of this source is comparable to 
that of GN20. The lack of submm emission could signal a 
warmer SED (see, e.g.. Chapman et al. 2008b). 

We notice that, due to the very conservative requirements 
applied, we finally retained only 3 galaxies of the 63 se- 
lected down to 30/iJy in GOODS-N with our radio/24/xm cri- 
terion as bona-fide z 3.5 starburst galaxy candidates. We 
believe that we efficiently rejected in this way all contami- 
nation from AGN dominated radio galaxies. Several of the 
discarded galaxies are though most likely to be starbursts at 
2.5 3.5, still very interesting objects but beyond the 

scope of the simple exercise described in this Section where 
we focus on the z > 3.5 redshift tail. We are likely also miss- 
ing several other genuine z ;^ 3.5 starburst, due to the very 
strict requirements that we have applied. E.g., that of a dis- 
tinct peak in IRAC beyond 7/im that can be weak or washed 
out due to the noise to the levels of our data; or that of X-ray 
non detections, powerful z 3.5 starbursts can well contain 
some unobscured AGN sometimes that doesn't necessarily af- 
fect much the radio emisson. Given the promising results of 
this attempt, we will defer to future work an investigation of 
the properties of the full sample of galaxies selected with this 
radio/24/im criterion. 

Although we by no means claim to have selected a com- 
plete sample, the existence of these three additional z > 3.5 
starburst candidates reinforces our conclusion that a substan- 
tial population of vigorous star-forming galaxies was already 
present at early epochs. 

7.2. Predictions for future Herschel surveys 

We now consider expectations for forthcoming Herschel 
surveys to detect very high redshift starbursts. Herschel 
will obtain FIR imaging with an unprecedented sensitivity at 
70/im, 100/im and with PACS (PogHtsch et al. 2006) 

and at 250/^m, 350/im and 52Q^im with SPIRE (Griffin et al. 
2007) . 

We computed the expected flux densities of GN20 in the 
Herschel bands, based on various assumptions. Using the 
best fitting z = 4.055 unsealed CEOl model shown in Fig. [8] 
we predict flux densities of about 40 — 50 mJy in the SPIRE 
bands, 24 mJy at 160/im, 4.8 mJy at lOO/im and 1.8 mJy at 
70/xm. The latter value is comparable, but still consistent, 
with the Spitzer 70/im upper limit. If these predictions are 
correct, galaxies similar to GN20 should be very easily de- 



tectable with Herschel, e.g., over the full 70 deg^ survey of 
the Hermes guaranteed time consortium'^, which could yield 
hundreds z > 4 hyperluminous infrared galaxies. Even if we 
used a colder template following Fig.[TOl normalized conser- 
vatively at a 1(7 lower 850/im flux density of 17 mJy, we 
would still predict flux densities of about 30 mJy at 350/im 
and 520/im, which is still above the expected flux limit of the 
70 deg^ Hermes survey. 

This suggests that, while strongly limited by confusion to 
quite bright flux densities (~ 15 mJy; Le Borgne et al. 2008), 
the Herschel SPIRE observations will efficiently detect a sig- 
nificant tail of very high redshift (z > 4) starbursts. Due to 
the relation between redshift and 850/^m flux density (Ivison 
et al. 2002; Pope et al 2006), Herschel might detect the bulk 
of the most active merger-powered galaxies in the distant uni- 
verse, and possibly accounting for an important fraction of the 
SFRD in the early universe. Clearly, a real challenge will be 
to recognize such very high redshift galaxies in the Herschel- 
selected samples, which will require deep multi-wavelength 
information. 

8. SUMMARY AND CONCLUSIONS 

The serendipitous detection of CO emission lines at z > 4 
from deep PdBl observations allowed us a first investigation 
of the nature and properties of vigorous star forming galax- 
ies at these redshifts. The main results of this paper can be 
summarized as follows: 

• CO emission lines are detected at the position of 
two bright SMGs in the GOODS-N field: GN20 and 
GN20.2a. A tentative CO signal is seen for an addi- 
tional radio galaxy in the field, named GN20.2b. 

• We identify the Unes of GN20 and GN20.2a as CO[4- 
3] at z = 4.05. A concordant identification is obtained 
when comparing optical photometric redshifts, ACS 
dropout color selection, radio-lR photometric redshifts, 
and optical spectroscopy from Keck. 

• Dust continuum at 3.3mm (0.65mm rest frame) is sig- 
nificantly detected for GN20, the first such a measure- 
ment in a distant starburst galaxy. 

• We show that reliable radio-IR photometric redshifts (to 
Az ^ 0.5) can be derived for these z = 4.05 galaxies us- 
ing only the photometry at 24/im, 850/im, 3.3mm and 
radio 20cm and CEOl models. This implies that local 
models can be used to provide a reasonable description 
of the z = 4.05 galaxies. 

• However, evidence is also found for these z 4 SEDs 
being somewhat colder than those of local galaxies of 
the same luminosities, confirming the results of P06 for 
lower redshift objects. In addition, the strong radio flux 
density of GN20.2a is taken as evidence for the pres- 
ence of an AGN in this galaxy. 

• The IR to radio SEDs suggest very high IR luminosi- 
ties of 2.9 and 1.6 x lO'-^Lg for GN20 and GN20.2, re- 
spectively. The CO emission lines imply molecular gas 
masses of 5 and 3 x lO'^Mp for GN20 and GN20.2, re- 
spectively. For GN20 we could derive an estimate of 
the dynamical mass of order of 2x10"^^^©, implying 

*^ |http : / / astronomy . Sussex . ac ■ uk/ jo/Hermes /| . 
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molecular gas fraction of order of 20% (or larger if the 
CO[4-3] transition is not thermalized). 

• The FIR to CO properties and SSFRs of these z > 4 
SMGs are similar to SMGs at lower redshift, 1 .5 < z < 
3. 

• Radio and submm stacking of B-band dropout galax- 
ies, together with the SED fitting of the ACS to IRAC 
SEDs, is used to show that typical z ~ 4 starforming 
galaxies have much lower SSFRs than SMGs at the 
same redshift. This suggests that the stellar mass-SFR 
correlation does not evolve much further in normaliza- 
tion beyond z ~ 2. 

• The GN20 and GN20.2 SMGs appear to live at the cen- 
ter of a strong sky overdensity of galaxies (with an ap- 
parent size of 1 — 2 comoving Mpc), and a very signifi- 
cant redshift spike is seen over the full GOODS-N area, 
consistent with a proto-cluster It is tantalizing to spec- 
ulate that this overdense environment might be playing 
a role in triggering such extreme star formation activity 
(see, e.g., Elbaz et al. 2007; Chapman et al. 2008b). 

• Applying the newly defined radio-IR photometric red- 
shift technique to the full sample of SMGs selected 
in the GOODS-N and Lockman Hole fields, we find 
that this provided estimates with a typical accuracy of 
Az/(1 +z) ~ 0.15. This is starting to approach the 
accuracy obtainable from optical photometric redshifts 
(Az/(1 +z) ^ 0.06; e.g., Brodwin et al. 2006). 

• Using this photometric redshift technique, we identify a 
small sample of galaxies in GOODS-N whose counter- 
part identification had likely failed in previous works. 
Lacking an obvious counterpart within the SCUBA 
beam, these SMGs could be at quite high redshift. 

• A total of up to eight SMGs in GOODS-N (includ- 
ing GN20, the two counterparts of GN20.2, GNIO, 
HDF850.1, GN13, GN28, GN04.2) are found consis- 
tent with being at z 4. This implies a space density 



of similar sources around 10^^ Mpc^^ and a sizeable 
contribution to the cosmic SFRD at z ~ 4 of order of 
lO^^M© yr"' Mpc"^. This is similar to that of UV- 
selected galaxies at the same redshifts. 

• We find that the GN20, GN20.2a and GN20.2b galax- 
ies and the most secure z > 4 SMG candidates satisfy 
the simple criterion Si 4 om. ^ SiA^un- By applying this 
criterion, with additional conservative constraints from 
optical, near-IR and mid-IR photometry, to a sample of 
VLA 20cm selected galaxies in GOODS-N, we iden- 
tify three additional galaxies that appear to be reliable 
z > 4 starburst galaxy candidates. A posteriori, two of 
the three are found to have a detection at submm/mm 
wavelengths as well. 

• Given the expected FIR fluxes of these z > 4 galaxies, 
future surveys with Herschel should allow us to charac- 
terize the properties of this population in detail. 
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